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ABSTRACT 

Fluctuations  in  the  horizontal  component  of  the  earth's 
magnetic  field  were  measured  at  a  land  site  (La  Mesa)  and 
on  the  floor  of  Monterey  Bay  at  a  depth  of  70  meters.   The 
measurements  span  a  19-hour  period  (August  17/18,  1982). 
Assuming  a  circular  polarization  basis  set,  analysis  was 
conducted  to  determine  the  coherence  between  the  two  sites 
in  the  frequency  range  .05-10  Hz.   Preliminary  analysis 
indicates  a  high  degree  of  coherence  from  .05-. 4  Hz.   Also, 
a  study  of  the  polarization  of  the  micropulsation  field  at 
each  site  was  conducted  in  terms  of  Stokes  parameters. 
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I.   INTRODUCTION 

This  thesis  research  is  part  of  an  ongoing  effort  being 
conducted  by  the  Naval  Postgraduate  School  to  obtain  a  long- 
terra  data  base  of  ELF  electromagnetic  noise  on  the  ocean 
floor.   The  result  of  this  project  should  be  a  better  under- 
standing and  interpretation  of  that  noise  and  its  sources. 

Magnetic  noise  is  important  to  the  geophysicist  and  to 
the  Navy.   Magnetic  techniques  are  applied  in  exploration 
for  minerals,  the  prediction  of  communications  propagation 
and  to  the  modeling  of  the  sphere  on  which  we  live.   Geo- 
magnetic studies  are  of  interest  to  the  Navy  for  their 
application  to  submarine  detection,  submarine  communications 
and  to  mine  warfare. 

The  objective  of  this  thesis  within  the  broader  context 
of  the  overall  project  was  to  design  and  validate  a  computer 
code  that  would  facilitate  data  manipulation  and  analysis. 
To  that  end,  the  present  computer  code  calculates  for  two 
sites  (one  on  the  ocean  floor,  the  other  on  land),  the 
individual  site's  power  spectral  densities  and  Stokes 
parameters  0,1,2  and  3  for  two  of  the  individual  site's 
orthogonal  components  in  the  horizontal  plane.   The  code 
also  computes  the  coherence  of  these  two  orthogonal  compo- 
nents and  their  phase  difference.   Lastly,  the  program 
computes  the  right  and  left  circular  polarization  power 
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spectral  densities  of  the  sites  and  the  coherence  between 
the  seafloor  and  land-based  sites. 

As  discussed,  the  present  code  applies  only  to  two 
sites  and  two  orthogonal  geomagnetic  components  at  each 
site.   However,  extension  to  several  sites  and/or  components 
should  be  an  easy  undertaking  due  to  the  style  of  the  present 
program. 
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II.   BACKGROUND 

A.   THE  MAIN  GEOMAGNETIC  FIELD 

Historical  observation  and  modern  measurement  demonstrate 
that  the  principal  source  of  the  geomagnetic  field  (approxi- 
mately 50,000  nanoteslas)  is  beneath  the  surface  of  the 
earth.   An  easy-to-picture  causal  interpretation  of  this 
phenomenon  is  painted  by  imagining  the  earth  with  a  dipole 
magnet  near  its  center  as  pictured  in  Figure  2.1.   The 
magnet  axis  is  tilted  at  an  angle  of  11.5  degrees  from  the 
earth's  rotation  axis.   This  model  provides  a  prediction 
accuracy  to  a  few  percent  of  the  measured  geomagnetic  field. 
Figures  accurate  to  less  than  one  percent  of  the  earth's 
measured  field  can  be  obtained  using  more  sophisticated 
models  [Ref .  1] . 

The  generally  accepted  physical  process  for  the  genera- 
tion of  the  geomagnetic  field  is  that  of  the  self-sustaining 
dynamo  [Ref.  2].   Figure  2.2  presents  a  simple  picture  of 
such  a  dynamo.   The  physical  requirements  of  this  theory 
are  mostly  met  in  the  interior  of  the  earth.   For  example, 
the  earth's  center  is  a  large  conducting  system.   Addition- 
ally, we  can  assume  that  the  earth's  molten  core  is  of  a 
differential  rotational  nature  with  turbulent  upwelling. 

Energy  sources  for  such  a  dynamo  are  most  likely  avail- 
able in  the  earth's  core  due  to  a  number  of  processes.   This 
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Figure  2.1    Eccentric  Dipole  Model  of  Geomagnetic  Field 
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Figure  2.2    Simple  Disk  Dynamo 

all  boils  down  to  mean  that  the  motion  of  a  conductor,  such 
as  the  earth's  molten  core,  in  a  magnetic  field  produces  a 
current  which  in  turn  produces  a  magnetic  field  in  support 
of  the  original  magnetic  field  [Ref .  3] . 

A  simplified  picture  of  the  dynamo  theory  originates 
by  assuming  that  in  the  distant  past  a  magnetic  "seed1*  field 
was  captured  by  the  earth's  conductive  core.   The  differen- 
tial rotation  of  this  core  "winds  up"  the  field  into  an 
intense  azimuthal  field,  such  as  would  happen  to  a  twisted 
rubber  band.   Turbulent  upwelling  of  the  core  tends  to 
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carry  this  azimuthal  field  outward,  converting  it  to  a 
toroidal  helical  field  within  an  associated  azimuthal 
circulating  current.   This  current  in  turn  generates  the 
earth's  dipolar  field. 

An  additional  advantage  of  the  self-sustaining  dynamo 
theory  lies  in  its  ability  to  explain  the  well -documented 
reversals  of  the  geomagnetic  field.   It  has  been  postulated 
that  two  or  more  coupled  "dynamos"  as  shown  in  Figure  2.3 
can  cause  reversals  of  the  main  geomagnetic  field  over 
geologic  time  scales  [Ref .  4] . 
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Figure  2.3    Twin  Disk  Dynamo 
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B.   GEOMAGNETIC  FLUCTUATIONS 

The  geomagnetic  field  changes  with  time.   The  changes  in 
the  field  have  periods  ranging  from  fractions  of  seconds  to 
thousands  of  years.   The  very  slow  variations  (periods  of 
years  or  more)  are  not  of  interest  here  but  are  known  as 
secular  variations  and  are  thought  to  be  geologic  in  origin. 

The  geomagnetic  field  changes  of  interest  here  are  those 
fluctuations  that  occur  over  the  period  of  a  day  or  less. 
It  is  customary  to  view  the  magnetic  field  as  the  stable 
main  field  with  these  comparatively  rapid  fluctuations 
superimposed  on  it.   Figure  2.4  presents  the  spectrum  of 
the  superimposed  fluctuations. 

The  fluctuations  (Figure  2.4)  that  are  superimposed  on 
the  main  field  range  in  intensity  from  .05  to  50  nanoteslas 
and  have  various  causes.   Some  major  contributors  to  these 
fluctuations  are  in  the  ionospheric  currents  overhead  and 
in  magnetic  storms  which  probably  originate  in  solar 
eruptions . 

Finally,  arriving  at  the  point  of  interest  here,  are 
those  disturbances  to  the  earth's  magnetic  field  that 
propagate  through  the  earth's  magnetosphere  as  hydrodynamic 
waves.   These  small  perturbations  cover  a  wide  range  in 
amplitude  and  duration  and  are  known  as  micropulsations . 
Figure  2.5  presents  a  graphical  display  of  these  disturbances 
Micropulsations  are  thought  to  be  the  result  of  fluctuations 
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Figure  2.4    Power  Spectrum  of  Geomagnetic  Disturbances 
Observed  on  the  Surface  of  the  Earth. 

[Claudis ,  Davidson,  and  Newkirk,  1971] 
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Figure  2.5    Field  Strength  of  Micropulsations 
[Campbell,  1966] 

in  the  incident  solar  wind  affecting  wave-particle  interaction 
in  the  magnetosphere  [Ref .  5] .   As  they  propagate  through 
the  conductive  layers  of  the  atmosphere,  they  are  trans- 
formed into  electromagnetic  waves  via  ionospheric  currents, 
which  we  detect  at  the  earth's  surface  as  magnetic  field 
fluctuations.   There  are  two  major  categories  of  micro- 
pulsations:   continuous  and  irregular. 

Smoothly  varying,  periodic  pulsations  are  classified  as 
(Pc)  micropulsations.   Their  distinct  periodicity  is  indica- 
tive of  a  resonance  inside  the  magnetospheric  cavity.   One 
model  proposed  to  explain  this  occurrence  of  smooth,  periodic 
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pulses  is  that  the  resonance  is  a  result  of  an  oscillatory 
hydromagnetic  wave.   This  wave  is  associated  with  groups 
of  ions  which  "bounce"  back  and  forth  in  the  magnetosphere 
[Ref.  6] . 

These  continuous  micropulsations  are  further  classified 
according  to  their  periods  as  follows: 

1.  Pel:   (0.15-5  second  period) 

Known  as  "Pearls,"  these  micropulsations 
are  generated  by  the  cyclotron  instability 
of  the  energetic  protons.   They  have  been  posi- 
tively correlated  with  solar  disturbances  and 
occur  during  daylight  hours  in  the  auroral 
zone  and  during  night  and  early  morning  hours  in 
the  midlatitudes .   Their  average  amplitude  is 
approximately  0.05-0.1  nanotesla  [Ref.  7], 

2.  Pc2:   (5-10  second  period) 

This  is  a  diurnal  phenomenon  that  shows  some 
positive  correlation  with  solar  activity  and  the 
seasons.   They  usually  decrease  in  their  period 
as  magnetic  activity  increases.   Their  average 
amplitude  is  0.1-1  nanotesla. 

3.  Pc3:   (10-45  second  period) 

This  micropulsation  is  similar  to  Pc2  pulsations 
except  for  the  difference  in  period. 

4.  Pc4:   (45-150  second  period) 

Sunspot  activity  appears  to  have  an  effect  on 
Pc4  pulsations.   The  Pc4  pulsations  also  change 
their  period  seasonally  with  an  average  amplitude 
of  5-10  nanotesla. 

5.  Pc5:   (150-600  second  period) 

These  large-scale  pulsations  typically  occur 
during  morning  and  evening  with  amplitudes  of 
10-100  nanotesla.   Their  duration  shows  a  strong 
geomagnetic  latitude  dependence. 

Irregular  (Pi)  micropulsations  appear  to  be  associated 

with  polar  magnetic  disturbances  and  are  observed  during 

night-time  hours.   The  theory  is  that  these  micropulsations 
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arise  from  disturbances  of  the  ionosphere  due  to  particle 
bombardment  from  the  magnetosphere.   Increased  solar  radia- 
tion levels  also  appear  to  have  an  effect  on  their  existence 
Irregular  micropulsations  are  classified  as  follows: 

1.  Pil:   (1-40  second  period) 

These  pulsations  have  a  small  amplitude  and 
usually  occur  at  night  and  early  morning. 
They  vary  in  intensity  from  0.01-0.1  nanotesla. 
They  demonstrate  a  positive  correlation  with 
auroral  disturbances. 

2.  Pi2:   (40-150  second  period) 

Pi2  pulsation  amplitude  range  is  on  the  order 
of  1-5  nanotesla.   They  usually  occur  during 
the  early  evening  hours  but  may  continue  through- 
out the  night.   The  periods  of  these  pulsations 
decreases  with  increasing  magnetic  activity. 

The  Pc5  and  Pc4  micropulsations  last  several  hours. 

However,  the  Pc3-1  micropulsations  and  irregular  pulsations 

have  a  maximum  duration  of  approximately  one  hour  but  may 

cease  after  only  a  few  minutes  [Ref.  8].   Above  3  Hz,  the 

most  pronounced  field  variations  are  the  earth  ionosphere 

cavity  resonances  of  Schumann  resonances  [Ref.  9] . 

C.   OCEAN  WAVE  EFFECTS 

In  addition  to  all  of  the  geomagnetic  fluctuations 
discussed  above,  if  working  in  the  ocean  environment,  one 
must  take  into  account  hydrodynamically  induced  magnetic 
fluctuations.   These  fluctuations  are  internally  produced 
resulting  from  the  movement  of  a  conductor  (seawater)  in 
an  external  magnetic  field.   The  effects  of  waves  have  been 
calculated  [Ref.  10]  and  are  pictured  in  Figure  2.6. 
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Figure  2.6 


Induced  Magnetic  Field  per  Meter  Amplitude 
of  the  Surface  Wave 

[Weaver,  1965] 
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Weaver's  equations  have  been  modified  to  incorporate  bottom 
reflection  and  transmission  with  the  result  that  the  effect 
of  bottom  decreased  the  magnitude  of  the  induced  magnetic 
field  by  less  than  3  dB  for  frequencies  greater  than  0.16  Hz 
in  water  depths  less  than  500  meters  [Ref.  11]. 
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III.   COLLECTION  SYSTEM 

A.   DATA  COLLECTION  SYSTEM 

The  data  acquisition  system  is  illustrated  in  Figure  3.1. 
Major  components  consist  of: 

(1)  coil  antenna  sensors  (2) 

(2)  preamplifiers  (2) 

(3)  signal  conditioners  (2) 

(4)  pulse  code  modulation  system  (1) 

(5)  radio  transmitter  (1) 

(6)  radio  receiver  (1) 

(7)  tape  recorder  (1) 

For  the  work  reported  in  this  paper,  two  systems  consisting 
of  components  (1)  through  (5)  were  deployed.   They  were 
similar  in  most  respects  with  the  major  exception  that  the 
ocean  equipment  was  enclosed  in  watertight  glass  spheres. 
A  thorough  description  of  the  data  acquisition  system  is 
given  in  Reference  12. 

1 .   Coil  Antenna  Sensors 

Each  coil  antenna  sensor  consists  of  two  continu- 
ously wound  coils  which  are  constructed  of  about  5460  turns 
of  18-gauge  copper  magnet  wire.   Glass  spheres  are  used  to 
maintain  the  watertight  integrity  of  the  coils  and  all  other 
immersed  equipment  for  the  ocean  configuration.   See  Figure 
3.2  for  the  ocean  system.   The  coils  are  mounted  orthogonally 
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Figure  3.2    Ocean  Data  Collection  System 
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in  a  nonmagnetic  frame  with  each  coil  having  a  resistance 
of  120  ohms  and  an  inductance  of  approximately  9.31  Henries. 

2 .  Preamplifiers 

The  preamplifiers  are  model  13-10A  low-noise  ELF 
amplifiers  manufactured  by  Dr.  Alan  Philips.   The  final 
stage  of  the  amplifier  contains  an  active  low-pass  filter 
which  provides  a  sharp  cutoff  for  frequencies  above  20  Hz. 
Frequencies  above  10  Hz  are  not  examined  in  this  paper.   The 
overall  amplifier  gain  for  inputs  of  less  than  2.5  milli- 
volts is  60  dB  in  power.   The  purpose  of  the  filter  is  to 
remove  60  Hz  noise  and  to  prefilter  the  data  before 
digitization. 

3.  Signal  Conditioners 

The  signal  conditioners  receive  analog  signals  from 
the  preamplifiers,  amplify  them  by  about  30  dB  and  limit 
signals  with  peak  amplitudes  of  7.5  volts  (modifiable)  from 
entering  the  pulse  code  modulation  system. 

4.  Pulse  Code  Modulation  System 

The  pulse  code  modulation  system  was  designed  and 
manufactured  by  Dr.  Robert  Lowe,  Lowecom,  Inc.   The  system 
offers  a  15-channel  analog  input  capability  with  selectable 
sampling  rates  of  2,  4,  8,  16,  32,  64  and  128  samples  per 
second.   The  PCM  system  uses  a  crystal  oscillator  and  inte- 
grated circuitry  to  develop  the  clocking  pulses,  a  16-channel 
analog  multiplexer,  a  16-channel,  12-bit  analog- to-digital 
converter  as  the  basis  for  the  pulse  coding.   The  crystal 
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clock  oscillator  loses  1  bit  in  1  million.   The  clock  pulse 
gates  the  analog  multiplexer,  the  analog-to-digital  converter 
and  additional  circuitry  that  form  the  pulse  code  words.   The 
PCM  output  is  organized  in  frames  with  each  frame  headed  by 
a  synchronization  word  which  is  followed  by  the  pulse-coded 
samples  from  PCM  channels  1  through  15. 
5.   Transmission  and  Recording 

After  encoding,  the  data  is  transmitted  via  VHF 
radio  to  a  Naval  Postgraduate  School  located  receiver  and 
recorded  on  an  instrumentation  tape  recorder.   For  the  data 
of  this  paper,  a  sample  rate  of  3  2  samples  per  second  was  used 

B.   DATA  ANALYSIS  EQUIPMENT 

1 .  Analysis  Hardware 

The  PCM  data  is  fed  into  a  PCM  decoder  which  places 
the  decoded  PCM  onto  a  9-track  800  BPI  computer  tape.   This 
computer  tape  is  read  into  the  Naval  Postgraduate  School  IBM 
3033  computer  for  manipulation  and  analysis. 

2 .  Data  Analysis  Software 

The  IBM  3033  is  used  to  convert  the  data  to  usable 
form  and  then  to  perform  any  operations  on  the  data  as  desired, 
At  present  the  computer  code  is  written  in  Fortran  IV  and 
converts  the  digital  data  to  voltage,  Fast  Fourier  trans- 
forms the  data,  calculates  power  spectral  densities,  deter- 
mines Stokes  parameters  0  through  3  and  provides  phase  and 
coherence  figures.   Appendix  A  contains  test  verification 
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of  the  code.   Appendix  B  contains  the  code  that  prepares 
data  from  one  site  for  storage  in  the  IBM  3033  Mass  Storage 
System  for  later  recall  and  comparison  with  another  site. 
Appendix  C  contains  the  code  that  prepares  and  manipulates 
the  code  from  a  second  site,  retrieves  from  the  Mass  Storage 
System  that  data  previously  prepared  and  manipulates  the 
two  data  sets  simultaneously. 
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IV.       COMPUTER   CODE 

A.  INTRODUCTION 

As  mentioned  earlier,  one  of  the  intentions  of  this 
effort  was  to  develop  a  code  that  would  permit  the  comparison 
and/or  simultaneous  geomagnetic  analysis  of  two  different 
sites.   The  present  data  collection  system  provides  for  the 
recording  of  transmitted  data  from  those  sites  onto  different 
tracks  of  the  instrumentation  tape  recorder.   Since  only  one 
track  of  that  analog  tape  can  be  placed  on  one  digital 
computer  tape,  a  means  of  storing  one  channel  had  to  be 
devised.   As  can  be  seen  by  a  quick  review  of  Appendices  B 
and  C,  the  present  technique  that  allows  this  simultaneous 
manipulation  of  data  from  two  sites  is  to  read  the  computer 
tape  from  one  site  and  to  process  that  data  to  the  point 
where  it  is  ready  for  analysis.   These  results  are  then  stored 
in  the  IBM  3033.   Subsequently,  the  same  routine  (with  modifi- 
cations to  fit  the  parameters  of  the  second  site)  are  applied 
to  data  read  from  the  second  site  computer  tape.   Now  both 
sets  of  data  are  in  the  same  form  and  ready  for  manipulation. 

B.  DATA  ANALYSIS  SOFTWARE 

The  two  programs  (one  for  each  of  two  sites)  are 
essentially  the  same  with  the  exception  that  one  program, 
program  B,  is  a  short  version  of  the  other,  program  A.   Thus, 
the  next  paragraphs  will  discuss  in  detail  the  longer  program. 
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This  description  fits  equally  well  to  program  B  except  that 
program  B  ends  after  applying  the  system  transfer  functions 
and  placing  these  results  in  the  Mass  Storage  System. 
1 .   The  Main  Program 

The  main  program  is  divided  into  sections  which 
perform  the  following  functions: 

*  Data  input 

*  Fourier  analysis  of  time  series  data 

*  Application  of  system  transfer  function.   Note:   At 
this  point,  program  B  reads  the  data  into  the  Mass 
Storage  System  and  ends.   Program  A,  at  this  point, 
retrieves  that  same  data  from  Mass  Storage. 

*  Data  averaging 

*  Calculation  and  plotting  of: 

(1)  power  spectral  densities  of  the  individual 
site  geomagnetic  orthogonal  components. 

(2)  orthogonal  component  coherence  of  the  geo- 
magnetic field  of  one  site. 

(3)  Stokes  parameters  0  through  3  of  the  individual 
site  orthogonal  components. 

(4)  single  site  right  and  left  circular  polarization 
power  spectral  densities. 

(5)  single  plane  coherence  of  right  and  left 
circular  polarization  between  sites. 

a.   Data  Input 

The  data  is  read  into  the  IBM  3033  via  the  sub- 
routine "read"  which  was  written  by  Tim  Stanton  of  the 
Naval  Postgraduate  School.   The  "read"  subroutine  reads  the 
PCM  data  from  the  computer  tape  and  converts  it  into  an 
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integer  from  0  to  4095.   The  main  program  than  converts 
this  integer  to  voltage.   See  Appendix  B  or  C. 

b.  Fourier  Analysis 

The  next  section  of  the  code  utilizes  a  sub- 
routine developed  at  the  Naval  Postgraduate  School  and 
called  "fourt."   "Fourt"  Fast  Fourier  transforms  the  voltage 
data  from  the  time  domain  to  the  frequency  domain.   See 
Appendix  A  for  validation  of  this  routine. 

c.  Transfer  Functions 

After  the  Fast  Fourier  transform,  the  program 
applies  the  transfer  functions  for  the  system.   Development 
of  the  transfer  functions  is  contained  in  Reference  12. 
The  transfer  functions  are  applied  up  to  10  Hz  as  straight 
line  functions.   The  data  is  now  in  nanoteslas. 

d.  Power  Spectral  Density 

This  quantity  is  computed  by  multiplying  the 
previously  calculated  spectral  density  summation  by  the 
observation  period  and  dividing  by  the  number  of  data 
groups  (typically  20)  that  the  computer  tape  was  segmented 
into,  providing  an  average  power. 

e.  Coherence 

The  next  steps  of  the  program  determined  for  one 
site  the  geomagnetic  orthogonal  components  Stokes  parameters 
0,  1,  2,  and  3,  and  the  phase  of  the  components  and  their 
coherence.   Next,  right  and  left  circular  polarization 
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parameters  are  computed  for  the  individual  sites  followed 
by  calculation  of  the  coherence  of  these  parameters  between 
the  two  sites. 

C.   PARAMETER  THEORY 

1 .   Stokes  Parameters 

A  useful  tool  for  determining  the  state  of  polariza- 
tion of  a  wave  are  the  Stokes  parameters  0,  1,  2,  and  3 
[Ref.  13]. 

Two  linearly  polarized  plane  waves  with  field 
vectors  in  the  directions  of  e.  and  e?   can  be  described  as: 


bl   £1  hl   e 


t     =  p  F  p  Ci  5-x  -  iwt) 
b2   £2  ^2  e 


These  waves  can  be  combined  to  give  a  general  plane  wave 
propagating  in  the  direction  k  =  kn : 

ft  (S,t)  -  (£j  Bx  *  £2  E2)  e^1  S>*  "  ^  (1) 


If  E,  and  E?  have  the  same  phase,  (1)  represents  a  linearly 
polarized  wave  with  its  polarization  vector  making  an  angle 

9  =  tan   (E  /E-)  with  E   and  with  a  magnitude  E  =  /E*  +  E*  , 
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If  E..  and  E  have  different  phases,  the  wave  (1)  is 
elliptically  polarized. 

In  the  simplest  case,  circular  polarization,  E,  and 
E2  have  the  same  amplitude  but  differ  in  phase  by  90°.  The 
wave  (1)  is  then: 

t    (J.t)  =  EQ  (£j  t  e2)  e  ***  -    iM«        (2) 

Then  the  components  of  (2)  are: 


Ex  (x,t)  =  Eq  cos  (kz  -  wt) 


Ey    (x,t)  =  EQ  sin  (kz  -  u>t) 


(3) 


Thus,  at  a  fixed  point  in  space,  the  fields  (3)  are  swept 
around  circularly  at  the  frequency  u>.   For  (e1  +is?),    the 
rotation  is  counterclockwise  if  the  observer  is  facing  the 
oncoming  wave.   This  wave  is  described  as  having  positive 
helicity.   (e,  -  ie-)  nas  negative  helicity  or  is  rotating 
clockwise  when  looking  into  the  wave. 

The  two  circularly  polarized  waves  (2)  form  an 
equally  acceptable  set  of  basis  fields  for  description  of 
a  general  state  of  polarization  by  introducing  the  complex 
orthogonal  unit  vectors : 
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e+  •  e,  =  0 
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Then  a  general  representation  of  (1)  is 


S  (S.t)  ■  (E+  S+  *  E.  J)  e  (i  k<X  "  iut)    (4) 


where  E+  and  E   are  complex  amplitudes.   If  E+  and  E  have 
different  amplitudes  but  the  same  phase,  (4)  is  an  elliptically 
polarized  wave  with  principal  axes  of  the  ellipse  in  the  plane 
defined  by  e,  and  e?. 

Thus,  the  polarization  content  of  a  plane  wave  is 
known  if  it  can  be  written  in  the  form  (1)  or  (4)  with  known 
coefficients  (E, ,  E_)  or  (E+ ,  E_) .   Conversely,  how  can  we 
determine  from  observations  on  a  wave  the  state  of  its 
polarization?   The  answer  lies  in  the  use  of  the  four  Stokes 
parameters.   Their  measurement  or  calculation  determines 
the  state  of  polarization  of  a  wave. 

For  a  wave  propagating  in  the  z  direction,  the  scalar 
products : 
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;2  .  a,  e2  .  g,  ej  .   £,  e*  -  2 


are  the  amplitudes  of  radiation  respectively,  with  linear 
polarization  in  the  x  direction,  linear  polarization  in  the 
y  direction,  positive  and  negative  helicity.   The  squares  of 
these  amplitudes  provide  a  measure  of  the  intensity  of  each 
type  of  polarization.   Phase  information  is  obtained  from 
cross  products.   If  we  assume  the  scalar  coefficients  of 
(1)  and  (4)  to  be  a  magnitude  times  a  phase  factor: 


„        16 1   „        i6 
El  =  al  e    >    E2  =  a2  e 


i6  +  „       ifi- 

E+  =  a_  e    ,  E_  =  a  e 


and  use  the  circular  polarization  basis  (e+,  e_) ,  we  have 


S   =  \i      .   E|2  +  \i      •  E|2  =  aj  +  a2 


S   =2Re[(e+  •  E  +  )   (£_  •  E)]  =  2a+  a_  cos(6_  -  6  +  ) 


S  =  2Im[(e+  •  E+)   (e_  •  EJ]  =  2a+  a_  sin(6_  -  5  +  ) 


S3  =  |e+  •  E|2  -  |e_  •  E|2  =  a2  -  a2 
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2 .   Coherence 

If  a  plane  wave  of  frequency  v  is  propagating  in  the 
z  direction,  it  has  components: 


ExCt)  =  a^t)^*'^  '  i2*ut 

that  are  mutually  orthogonal  and  at  right  angles  to  the 
direction  of  propagation  [Ref .  13] .   Then  the  complex 
degree  of  coherence  can  be  written  as : 


J 


xy 


xy   /D —  /J — 
xx   yy 


where   J      ....    are    the   elements   of   the   coherency  matrix 

xx' 


3  = 


E   E  E   E 

xx  x  y 


y  x 


* 

E   E 

y  y 


aiaze 


i<f>i  -  i  4>  2 


aia2e  a2 


The    complex   factor   y        is    the    complex   degree   cf   coherence   and 

xy 

is  a  measure  of  the  correlation  between  components  of  the 
field  vector  in  the  x  and  y  directions. 
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D.   SOFTWARE  TESTING 

The  basis  set  for  the  code  calculations  is  that  of  a 
circularly  polarized  wave.   Thus  a  circularly  polarized 
wave  was  chosen  as  the  test  input.   The  wave  was  created  by 
generating  a  4  Hz  sine  wave  for  sites  A  and  B  x  coil  input 
and  4  Hz  cosine  wave  for  sites  A  and  B  y  coil  input.   The 
results  of  this  test  are  valid  only  at  the  input  frequency 
as  a  measure  of  the  code  veracity.   One  of  the  characteris- 
tics of  any  FFT  routine  is  that  small  residuals  are  created 
throughout  the  spectrum  of  interest.   These  residuals  are 
primarily  caused  by  the  finite  observation  time.   These 
small-valued  artifacts  will  be  carried  throughout  all  subse- 
quent parameter  calculations  and  will  appear  in  the  test 
plots.   This  inclusion  of  unexpected  values  at  frequencies 
other  than  the  test  input  of  4  Hz  should  be  considered  as 
the  result  of  the  artificiality  of  the  single  value  test 
input.   Input  of  real  world  data  that  covers  a  full  spectrum 
of  values  obscures  the  "FOURT"  routine  propensity  to  create 
these  artifacts. 

1 .   Sites  A  and  B  Orthogonal  Component 
Power  Spectral  Density  Test 

Spectral  density  is  defined  as: 

S   =  i  E  (x.)  (x.) 
x   N  ,    l    l 

for  site  A  x  coil  where  x  is  the  Fourier  Transform  of  some 
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input  to  the  x  coil.   It  follows  that  the  average  power 
spectral  density  is: 


PSD 


xa 


Time 

N 


I    (XaHXa)*  =  Sxa  ■  Time 


where  time  is  the  observation  time.   Thus,  for  the  test 
inputs  of  sine  4  Hz  into  the  x  coils  at  sites  A  and  B  and 
the  cosine  4  Hz  into  the  y  coils,  we  expect  and  obtain  (see 
Figures  A.l,  A. 2,  A. 11,  A. 12)  a  peak  magnitude  of  .5  at 
4  Hz  at  each  site. 

2 .   Sites  A  and  B  Orthogonal  Component  Coherence  Test 
Coherence  as  defined  above  is : 
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for  two  orthogonal  components  where  S    is  defined  by 

o  r  Xya 


xya 


TT1*  0<a)(Ya)* 


With  X  and  Y  our  transformed  inputs  and  considering  the  same 
inputs,  we  expect  to  see  a  coherence  of  magnitude  1  at  a 
frequency  of  4  Hz  and  its  harmonics.   Figures  A. 8  and  A. 18 
demonstrate  the  validity  of  the  coherence  calculations  of 
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the  code.   Also  the  value  of  1  for  the  test  case  coherence 
at  frequencies  other  than  4  Hz  is  a  manifestation  of  the 
residual  effect  discussed  above. 

3.  Sites  A  and  B  Stokes  Parameter  1  Test 
Using  the  definitions  developed  in  preceding 

paragraphs,  Stokes  1  is: 

Sla  =  2a+a_sin  (6_  -  6 +)  =  2 [Sx  -  S  ]  =  0 

for  S   =  S  as  is  the  case  for  this  test  input  and  as  is 
x    y 

shown  in  Figures  A. 4  and  A. 13. 

4.  Sites  A  and  B  Stokes  Parameter  2  Test 

Again  using  definitions  developed  in  the  preceding 
paragraphs,  Stokes  parameter  2  is: 

soo  "  2a  a  sin  (6   -  6X)  =  0  for  (6   -  6 J  =  0 

Figures  A. 5  and  A. 15  validate  the  code  for  this  calculation 

5 .  Sites  A  and  B  Stokes  Parameter  3  Test 
Stokes  3  is  : 


S3a  =  a+  "  a-  =  -(-707)2  -C707)2  =  -1 


Figures  A. 6  and  A. 16  validate  the  code  for  this  calculation 
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6.   Sites  A  and  B  Stokes  Parameter  0  Test 
Stokes  0  is  : 


Srt  =  a?  +  a2  =  2  [S   =  S  . ]  =  2 
oa    +    -      L  xa    xbJ 

(See  Figures  A. 3  and  A. 13.) 

7.  Sites  A  and  B  Phase  Difference  of  x  and  y  Test 

The  phase  difference  of  sine  4  Hz  versus  cosine  4  Hz 
is  90°  or  1.57  radians.   (See  Figures  A. 7  and  A. 17.) 

8 .  Sites  A  and  B  Two  Orthogonal  Component  Power 
Spectral  Density  with  Right  and  Lett  Circular 
Polarization  Basis  Set  Test" 

The  development  of  the  spectral  densities  for  this 

case  is  similar  to  the  development  of  spectral  density  for 

a  single  component.   Starting  with  the  cross  spectra  density 

definition  for  right  and  left  circular  polarization 


BP    =  (X  +  iY)     (right  polarized) 


BM    =  (X  -  iY)     (left  polarized) 
xya    i.      -»     v     r 


It  follows  that  the  spectral  density  is: 


SP    =  i  Z  (BP   ) (BP   )      (right  polarized) 
xya   N  1  v  xya' ^  xya' 


-i  n  * 

SM    =  ir   Z  (BM    )  (BM    )      (left  polarized) 
xya   N  ,  v  xya'  *•  xya' 


43 


The  power  spectral  densities  will  be: 


SP    =  ±±2E.  i   SP 

xya    N   -,   xya 


SM    =  li^-  I  SM 
xya    N   ,    xya 


where  the  time  is  the  observation  time.   Examples  of  this 

calculation  with  a  test  input  of  4  Hz  right  circularly 

polarized  wave  are  in  Figures  A. 9,  A. 10,  A. 19  and  A. 20. 

9 .   Coherence  of  Sites  A  and  B  Right  and  Left  Circular 
Polarization  Basis  Set  Test" 

Right  and  left  circular  polarizations  are: 


SP  . 

COP  =  aoxy_     (right  polarization) 

/SP~  /3MT 
a    b 


SM 

COM  =  — — —    (left  polarization) 

/SM~  /5MT 
a    b 


where  SP  .    and  SM  ,    are  the  cross  spectra  densities  of 
abxy       abxy 

the  two  orthogonal  components  at  the  sites  A  and  B. 


T-    n 
SP  ,    =  lime  z  (BP   )  (BP    )     (right  polarization) 
abxy    N   n  L  xyaJ  L  xyaJ     ^   6   F 
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SMabxy  =  TT*  I    CBMxya)  CBMxyb^     <le£t  PO^rization) 

The  test  case  of  sine  4  Hz  as  the  xa  and  xb  inputs  versus 
cosine  4  Hz  as  the  ya  and  yb  inputs  yields  a  coherence  of  0 
for  left  circular  polarization  and  1  for  right  circular 
polarization.   (See  Figures  A. 20  and  A. 21.) 
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V.   ANALYSIS 

Preliminary,  "eyeball"  analysis  was  conducted  on  data 
collected  over  a  period  from  1122  local  17  Aug  82  to  0639 
local  18  Aug  82.   The  timeframes  inspected  within  this 
window  were : 

TIME  (local)  DATE 

1122-1254  17  Aug  82 

1301-1408  17  Aug  82 

2028-2200  17  Aug  82 

0121-0251  18  Aug  82 

0259-0431  18  Aug  82 

0507-0639  18  Aug  82 

Original  intentions  were  to  collect  data  over  a  full  24-hour 
period.   However,  the  batteries  at  the  Monterey  Bay  site  went 
dead  after  20  hours.   Of  the  20  hours  of  collected  data, 
eight  hours  appeared  dominated  by  system  noise  and  were  not 
included  in  the  analysis.   However,  the  remaining  data  of 
approximately  eight  hours  spans  a  24-hour  period  and  provides 
evidence  of  some  obvious  trends  that  lend  themselves  to  the 
cursory  analysis  conducted  here.   Representative  plots  of 
the  calculated  parameters  are  in  Appendix  D.   Technical 
Report  NPS  61-83-005  contains  a  complete  set  of  plots  of  the 
calculations  that  span  the  full  20  hours  of  data  collection. 
This  complete  set  of  plots  was  used  in  the  analysis  that  follows 
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A.  POWER  SPECTRAL  DENSITIES 

Both  sites  exhibit  a  marked  decrease  in  power  spectral 
density  with  frequency.   There  is  an  obvious  exception  to 
this  observation  at  4  Hz  at  the  land  site.   The  sharp  peak 
at  4  Hz  is  attributed  to  60  Hz  aliasing  which  is  a  result 
of  the  Fast  Fourier  Transform  (FFT)  subroutine.   Fast 
Fourier  algorithm  construction  is  such  that  data  being 
transformed  may  be  convolved  back  onto  itself  [Ref.  15]. 
This  convolution  is  routinely  avoided  by  using  a  sample 
rate  that  is  twice  the  upper  frequency  of  the  domain  of 
interest.   A  sample  rate  of  32  samples  per  second  was  used 
in  this  experiment  and  should  have  precluded  aliasing  of 
the  frequency  of  60  Hz,  especially  when  combined  with  the 
amplifier  filter  cutoff  of  20  Hz.   However,  the  60  Hz 
component  of  the  magnetic  field  was  measured  and  found  to 
be  4  nanotesla/Hz  which  is  105  higher  than  magnetic  noise 
data  [Ref.  16].   The  20  Hz  filter  cannot  completely  filter 
out  signals  of  this  magnitude.   Thus  the  4  Hz  signal  is  a 
result  of  aliasing  of  the  strong  ambient  60  Hz  signal. 

B.  STOKES  PARAMETERS 

All  Stokes  parameters  demonstrate  reasonable  agreement. 
That  is,  Stokes  0  is  obviously  a2.  +  a2  and  Stokes  3  is 
obviously  a2.  -  a2  as  shown  in  Figures  D.3,  D.6,  D.13  and  D.16 
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1 .  Ocean  Site 

The  ocean  site  Stokes  2  and  3  parameters  exhibit  a 
time  dependence.   That  is,  in  the  early  a.m.  hours  the  large 
magnitude  of  Stokes  2  is  the  noticeable  feature.   This 
magnitude  decreases  at  sunrise  and  through  the  day  until  it 
reasserts  itself  after  sunset.   A  review  of  Stokes  parameters 
2  and  3  over  the  total  observation  time  (see  Technical  Report 
NPS-61-83-005)  indicates  that  this  falloff  is  independent  of 
frequency.   The  increase  in  magnitude  of  Stokes  2  after 
sunset  is  also  independent  of  frequency  with  the  exception 
of  some  very  discrete  spectral  components. 

2 .  Land  Site 

The  land  site  Stokes  2  and  3  parameters  show  a  time 
dependence  with  a  shorter  cycle.   That  is,  the  magnitude  of 
Stokes  2  is  large  near  sunset,  followed  by  its  decrease 
during  the  morning  hours.   The  early  afternoon  hours  reflect 
an  increase  in  Stokes  2  until  after  sunset.   Figures  5.1 
through  5.4  show  the  decreased  magnitude  of  Stokes  parameters 
2  and  3  for  the  morning  hours. 

C.   X  AND  Y  COIL  PHASE  DIFFERENCE 
1.   Ocean  Site 

The  use  of  Stokes  parameters  in  the  calculation  of 
the  phase  difference  dictates  that  changes  in  the  time 
series  averaged  Stokes  parameters  should  be  evident  in  the 
x  and  y  component  phase  difference.   This  is  the  case  as 
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shown  in  Figures  D.7  and  D.17.   The  phase  difference  from 
sunset  to  sunrise  shows  a  definite  decrease  with  an  "eye- 
ball" average  of  .2  to  .4  radians.   Daylight  hours  show 
a  difference  ranging  from  .5  to  .3  radians  with  the  greatest 
difference  near  local  noon. 
2.   Land  Site 

The  land  site  phase  difference  reflects  its  Stokes 
parameters  with  a  difference  of  about  .3  radians  in  the 
higher  frequencies  at  all  hours  and  a  difference  of  .1 
radians  during  the  daylight  hours. 

D.   X  AND  Y  COIL  COHERENCE 

1 .  Ocean  Site 

Time  dependence  is  again  a  strong  factor  with  high 
degrees  of  coherence  in  the  ocean  site  x  and  y  coils  being 
greatest  during  hours  of  darkness  and  diminishing  at  sunrise. 
Figures  5.5  and  5.6  are  examples  of  this  phenomenon.   The 
exception  to  this  observation  is  at  local  noon  when  the 
coherence  of  the  ocean  coils  is  again  very  high  but  diminishes 
before  and  after  the  noon  time  data. 

2 .  Land  Site 

The  land  site  x  and  y  coil  coherence  reflects  the 
previous  parameter  cycle  time  dependence.   There  is  a  high 
degree  of  coherence  in  the  lower  frequencies  at  sunrise 
that  is  followed  by  diminishing  coherence  in  the  higher 
frequencies.   Local  noon  shows  an  improved  coherence  at 
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these  frequencies  but  is  followed  by  a  degraded  coherence 
after  sunset. 

E.  CIRCULAR  POLARIZATION  POWER  SPECTRAL  DENSITY 
There  is  no  marked  contrast  between  right  and  left 

circular  polarization  power  spectral  densities.   Both 
calculations  provide  similar  results. 

F.  COHERENCE  OF  CIRCULAR  POLARIZATION 

No  preference  for  right  versus  left  circular  polarization 
is  discernable .   The  most  obvious  trend  is  that  of  a  high 
degree  of  coherence  in  the  frequency  range  .05  Hz  to  .4  Hz 
as  shown  in  Figures  5.7  and  5.8.   Generally,  coherence  of 
the  two  sites  decreases  thereafter  with  frequency. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

In  view  of  the  limited  data  base  and  the  preliminary 
state  of  the  analysis,  only  tentative  conclusions  are 
possible  at  this  point. 

(1)  In  the  frequency  range  .05  to  .4  Hz,  the  ocean  and 
land  data  exhibit  a  high  degree  of  coherence.   At  higher 
frequencies  the  coherence  rapidly  diminishes. 

(2)  Both  polarized  and  unpolarized  signals  were  observed 
with  a  preponderance  of  high  degrees  of  polarization  occurring 
at  the  lower  frequencies.   The  nature  of  the  polarization 
appears  to  be  a  function  of  both  time  and  frequency. 

B.  RECOMMENDATIONS 

The  first  and  obvious  recommendation  is  to  collect  more 
data.   This  is  already  in  the  planning  stage.   The  second  and 
also  obvious  recommendation  is  to  conduct  a  more  thorough  ana- 
lysis of  the  data  on  hand.   This  also  is  already  being  planned 
More  concrete  recommendations  concern  data  processing. 

One  90-minute  real  time  data  tape  requires  12  computer 
cylinders  of  storage  space.   Thus,  for  the  approximately 
eight  hours  of  data  computed  and  analyzed  here,  72  cylinders 
of  storage  were  used.   This  is  too  much  and  is  mostly  a 
matter  of  the  computer  code  employed.   A  little  work  can 
provide  a  program  that  would  require  less  computer  storage 
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space.   This  will  be  a  real  necessity  as  a  data  base  of 
weeks  or  months  is  compiled  for  analysis.   It  is  recommended 
that  the  present  code  by  "streamlined"  to  use  the  IBM 
3033  as  a  calculation  "scratch  pad."  Then  a  few  lines  of 
"Computer  Job  Control  Language"  can  be  employed  to  store 
the  calculations  on  a  high-density  tape  for  future  use  at 
the  Naval  Postgraduate  School  or  anywhere  in  the  world. 
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APPENDIX    A 

TEST     PROGRAM 
//JOEM1MB    JOB     (1592,01  29) ,' FISHER    SMC.  1 399 • , CLASS  =  C 

//♦MAIN    ORG=NPGVM1.  159  2P, LINES= ( 60) 

//♦FORMAT    PR,DDNAME=PLOT.SYSVECrR,DEST=LOCAL 

//    EXEC    FRTXCLGP,PARM.  LKED=' LIST , MAP, XREF' , REGION. GO=2700K 

//SYSUT1  DD       UNIT=SYSDA,SPACE=(CYL,  (8,8)  ) 

//SYSUT2  DD       UNIT=SYSDA,SPACE=  (CYL,  (8,8)  ) 

//SYSLIN  DD       SPACE= (6080  ,  (80,30)  )  ,UNIT  =  SYSDA 

//FORT.SYSIN   DD    DSN=MS S. SYS3 . NON IMSL. SOURCE (FODRT) , DIS?= SHR 

//  DD       * 

C  A   LARGE    NUMBER    OF    ARRAYS    ARE    UTILIZED 

C  WITH    THE    INTENTION    OF    MAINTAINING    EACH 

C  BIT    OF    ORIGINAL    DATA    OR    CALCULATION, 

C  BEADY    FOR    RECALL    AT    &NY    TIME    FOR    FURTHER 

C  COMPUTAITON.    EQUIVALENT    ARRAYS    ARE    EMPLOYED 

C  BECAUSE   THE    PLOTTING    ROUTINE    'DRAWP'     IS 
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NOT    ABLE    TO    HANDLE    COMPLEX    ARRAYS. 
ARRAYS    ■  HB1,  'RTB',  '  &LA3«     AND    TITLE'     ARE 
USED    TO    GENERATE   THE    VERSATEC    PLOTTER    OUTPUT, 
OCEAN    ARRAYS    ARE    INDICATED    VIA    THE    SUFFIX 
'0',     AND    LAND    ARRAYS    ARE    INDICATED    WITH    'L'. 
C0MPLEX*8    XX0(3  192)  ,YYO(3192)  , 

C    CSO  (8192)  ,C00(8  192)  ,  ZXD  (3192)  , 

C    CRO  (8192)  ,CT0(3  192)  ,CUD  (3192)  , 

C    CVO  (819  2)  ,CW0(3  19  2)  fCPO  (319  2)  , 

C    ZYO (819  2)  ,C0M3L  (819  2)  ,C3L  (8192)  , 

C    CUL  (819  2)  fCVL(3  192)  rCWL(8192)  , 

C    COPOL(8192)  ,  XXL  (819  2)  ,YYL  (3  192)  , ZXL(819  2)  , 

C    ZYL  (8192)  ,CSL(3  192)  ,BPD  (3192)  , 

C    BMO  (819  2)  ,3PL(8  192)  ,  3ML  (8192)  , 

C    SPOL  (8192)  ,SM3L  (819  2)  ,SPD  (8  192)  , 

C    SPL  (8192)  ,SM0(8  192)  ,  CTL  (3192)  , 
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C   CPL  (8192)  ,CRL(3  192)  ,SML(3192) 


DIMENSION    TI  MB  (8192)  ,  FRE2  (8  192)  , 


C    WORK  (16384)  ,FRQ2(81  92) 


INTEGER*4    IT  B  ( 1  2) /1  2*0/ 


REAL*4  RTB(28)  /28*0.0/, 


C    CROO  (16  3  84)  ,CFD 0(16 3  84)  , 3U0O  (16  384)  , 


C    COOO  (16  3  8  4)  ,  CrDO(  16  3  84)  ,  3V0  0  (16384)  , 


C    CWOO  (16384)  ,CTLL(1S384)  , 


C    CPLL  (16  3  84)  ,  CRLL(  16  3  84)  ,  3DLL  (16384)  , 


C    CULL  (16384)  ,C7LL(16384)  , 


C   CHLL  (16384)  r  COM  OLL  ( 16334)  ,  COPOLL  (16384)  , 


C    SPOO  (16384)  ,SPLL(16384|  , 


C    SMOO  (16384)  ,SMLL(15384) 


REAL    ALAB(22)  /•  PDXO'  ^PDYO^'STIO1, 


C    'ST20' , ■ ST30', ' ST0O» , 'PHAO' , 


C    'GAMC  ,  ' PDXL',  '  PDYL* , ' ST1 L '  , • ST2L » , 
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C    « ST3L' , ■ STOL',  ■  PHAL* , '3AML' , 

C    '  SPLO1 , • SPLL1, • SHOO1 ,  'SMLL1 , f CPOL», »CMOL'/ 

REAL*8    TITLE  (12) 

EQUIVALENCE  (TITLE  (1)  ,RTB(5)  )  ,  (CROO(1)  , 
C    CEO  (1)  )  ,  (CP30(1)  ,  CPO  (1>  )  , 

C  (COOO(1)  ,COO  (1}  )  ,  (CTOO(1|  ,  CTO(1)  )  ,  (CUOO  (1)  ,CUO  (1)  )  f 
C  (CWOO(1)  ,CWO  (1)  )  ,  (CVOO(1)  ,CVO(1)  )  ,  (CTLL(1)  ,CTL(1)  )  , 
C  (CPLL(1)  ,CPL(1)  )  ,  (CRLL(1|  ,CBL(1)  )  f  (COLL(1)  ,COL(1)  )  , 
C  (CULL(1)  ,CUL  (1)  )  ,  (CVLL(1|  ,CVL(1)  )  ,  (CWLL(1)  ,CWL  (1)  )  , 
C  (COPOLL  (1)  ,COP3L(1)  )  ,  (COMOLL(1)  , 
C  COMOL(1)  )  ,  (SPD3  (1)  ,SPO(1)  )  , 
C     (SPLL(1)  ,  SPL  (1|),  (SMOO(1|  ,  SMO  (  1)  )  ,  (SMLL  ( 1)  ,SML  (1)  ) 

DATA    XXO ,YYO,CSO,COO,ZXO,ZYO,CRO,CTO,CUO,CVO,CWO, 
C    XXL, YYL,ZXL,Z*L  ,CSL, BPO, BMO ,BPL, 
C    BML,  SPOL,SMOL,SPO,SPL,SMO, 
C    SML, CTL,CRL,CDL,CVL,COPOL,CWL, 
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C   COPOL,COMOL/278528*  (0.0,0.0)/ 

DATA   TIME,FREQ, FRQ2/24576* 0. 0/ 

IFRAME=8192 

NR=30 

FNR  =  FLOAT  (NR) 

P=0.  0 

DO    7  0   L1=1,NR 

THE    DO   L03P    ENDING    WITH    STATEMENT    70    ENABLES 
THE    PR03RAM    TO    PROCESS     A    LARGE    AMOUNT    OF    DATA 
EY    REPEATING    THE    PROCESS    IN    BLOCKS. THE    DATA 
POINTS    FROM     EACH    RUN    THROUGH    THE    DO    LOOP    ARE 
ADDED   TOGETHER    AND    EVENTUALLY    AVERAGED    BY    THE 
NUMBER    OF    RUNS    THR0U3H    THE    DO    LOOP.     '  NR« 
REPRESENTS     THE    NUMBER    OF    DATA    SEQUENCES    TO 
BE    AVERAGED.     1    SEQUENCE    CURRENTLY     EQUALS 
8192    DATA    POINTS    FDR    EACH    CHANNEL    OR    256 
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SECONDS    OF     DATA.     FOR    THIS    TEST    SINE    AND 
COSINE    WAVES    OF    U HZ    WIL    BE    USED. 

DO    60  JJ=1,IFRAME 

XXO  (JJ)  =  SIN(P*.  78  53  98163) 

YYO  (  JJ)  =COS(P*.  785398163) 

XXL (JJ)  =SIN(P*.  785398163) 

YYL  ( JJ)  =COS  (P*.  785398163) 

P=P  +  1 
60   CONTINUE 

THE    FOLLOWING    SECTION    GENERATES    THE    TIME 
AND    FREQUENCY    ARRAYS    AND    NORMALIZES 
THE    INPUT    PCM    DATA    TO    VOLTAGE    FORM    IN 
PREPARATION     FOR    THE    FAST    FOURIER    TRANSFORM 
TO    THE    FREQUENCY    DOMAIN. 

N=8  1  92 

FN=FLOAT  (N) 
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DELTAT= 1./32. 

T=FN*DELTAT 

DELTAF=1./T 

DO    2C   J  =  1,N 

TIME  (J)  =DELT  AT*  FLOAT  (J| 

FREQ  (J)  =DELTAF*FLOAT  (J) 

FRQ2  (J)  =ALOG10  (FREQ  (J)  ) 
20      CONTINUE 

THE     NEXT    FOUR    STATEMENTS    PERFORM 
AN    FFT    ON    THE    INPUT    TIME    SERIES 
DATA. 

CALL    FOURT(XXD, N, 1, -1, 0,tfORK) 

CALL    FOURT(YYO, N, 1, -1,0, WORK) 

CALL    FOORT(XXL,N, 1,-1,0, WORK) 

CALL    FOORT(YYL,  N, 1,-1,0,  WORK) 

DO    UO   K4=1,N 
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XXO  (K4)  =XXO(KU)  /FN 
YYO  (K4)  =YYO(K'4|  /FN 
XXL (KU) =XXL(K4>  /FN 
YYL  (K4)  =YYL(K4)  /FN 
40      CONTINUE 


THE    NEXT    LOOP    CALCULATES    THE 


UNNOFMALIZED    SPECTRAL    DENSITIES 


FOR     SINGLE     SITE    ORTHOGONAL 


COMPONENTS     OF    THE    GEOMAGNETIC    FIELD, 


THE    INDIVID  UAL    SITS    OROS    SPECTRA 


EETWEEN    COMPONENTS,     THE    INDIVIDUAL 


CRTHOGONAL    COMPONENTS    OF    RIGHT    AND 


IEFT    CIRC3L  AP.LY    POLARIZED     SPECTRA, 


AND    TWO    SITE    CROSS    SPECTRA    FOR    RIGHT 


AND    LEFT    CIRCULAR    POLARIZATION, 


DO    3C   11=1, N 
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ZXO  (II 


Z  YO  ( II 


CSO  ( II 


ZXL (II 


Z  YL  ( II 


CSL (II 


BPO  (II 


BMO  (H 


BPL (II 


BML (II 


=  ZXO(II)  +  (XXO  (II)  *C3NJG  (XXO  (II)  ) 


=  ZYO(II)  +  (YYO(II)  *CONJG  (YYO(II)  ) 


=  CSO(II)  +  (XXO  (II)  *C0NJG  (YYO  (II)  ) 


=  ZXL (II)  ♦  (XXL  (II)  *CDNJG  (XXL  (II)  ) 


=  ZYL(II)  +  (YYL(II)  *CONJG  (YYL(II)  ) 


=  CSL  (II)  +  (XXL  (II)  *CTDNJG  (YYL  (II)  ) 


=  XXO  (II)  +(  (0.  ,1.)  *YYO  (II)  ) 


=  XXO(II)  -(  (0.  ,1.)  *YYO  (II)  ) 


=  XXL (II)  +  ( (0. , 1.)  * YYL  (II)) 


=  XXL(II)  -  (  (0.  ,1.)  *YYL  (II)  ) 


SPOL  (II)  =SPOL(I  I)  ♦  (BPO  (II)  *CONJG  (BPL  (II)  )  ) 


SMOL  (II)  =  SMOL(II)  +  (BMO  (II)  *CONJG  (BML  (II)  )  ) 


SPO  (H)  =SPO(II)  +  (BPO  (II)  *CONJG  (BPO  (II)  )  ) 


SPL  (II)  =SPL  (II)  +  (BPL(II)  *CONJG  (BPL(II)  )  ) 


SMO  (H)  =SMO(II)  +  (BMO  (II)  *CONJG  (BMO  (II)  )  ) 


SML  (II)  =SML(II)  +  (BML  (II)  *CONJG  (BML  (II)  )  ) 
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3  0      CONTINUE 
70      CONTINUE 


THE    NEXT    LOOP    NORMALIZES    THE       ABOVE 


SPECTRA    AND    CALCULATES    POWER    SPECTRAL 


DENSITIES. 


DO    33   13  =  1, N 


ZXO  (13)  =ZX0(I3)  *T/FNR 


ZYO  (13)  =ZY0(I3)  *T/FNR 


CSO  ( 13)  =CSO  (13)  *T/FNR 


ZXL  (13)  =ZXL(I3)  *T/ENR 


ZYL  (13)  =ZYL(I3)  *T/FNR 


CSL  (13)  =CSL(I3)  *T/FNR 


SPOL  (13)  =SPOL(I  3)  *I/?NR 


SPO  (13)  =SP0(I3)  *T/FNR 


SPL  (13)  =SPL(I3)  *T/FNR 


SMO  (13)  =  SM0(I3|  *T/FNR 
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SML  (13)  =SML(I3)  *T/FNR 
SMOL  (13)  =SMOL(I  3)  *T/FNR 
33      CONTINUE 

THE    NEXT    LOOP    CALCULATES    STOKES 
PARAMETERS    0    THROUGH    3    OF    THE 
INDIVIDUAL     SITE    ORTHOGONAL 
COMPONENTS     AND    THE    COHERENCE    OF 
THE    PLANAR     CIRCULAR    POLARIZATION 
PARAMETERS     BETWEEN    SITES. 
DO    44   14  =  1, N 

CTO  (14)  =  (ZXO  (I '4)  +  ZYO(I4)  |  *2./T 
CPO  (14)  =  (ZXO  (I'4)-ZYO  (14)  )  *2./(T*CTO  (14)  ) 
CRO  (14)  =  (4*CS0  (14)  )  /(CIO(I4)  *T) 

COO  (14)  =CSO(I4)  /(CSQRT(ZXO (14)  )  *CSQRT  (ZYO  (14)  ) ) 
CWO  (14)  =  ATAN2(AIMAG  (COO  (I  4 )  )  ,  RE  AL  (COO  ( 14)  )  ) 
COO  (14)  =CSQBT(COO  (14)  *CDS JG  (COO  (14)  )  ) 
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CTO  ( 14 


COO (14 


C  VO  ( 14 


CTL  ( 14 


CPL (14 


CRL  (14 


COL  ( 14 


CTL  ( 14 


CUL  (14 


CVL  ( 14 


CWL (14 


=4.34  29  4 48 *C LOG (CTO  (14) ) 


=4.3429448*CLOG(ZXD  (14)  ) 


=  4.34  29  4 48 *CLOG(ZYO  (14)  ) 


=  (ZXL  (14  )  +ZYL  (14)  )  *2./T 


=  (ZXL  (14)  -ZYL  (14}  I  *2./(T*CTL  (14)  ) 


=  (4*CSL (14) ) /(CTL (14)  *T) 


=CSL  (14)  /(CSQRT(ZXL  (14)  )  *CSQRT  (ZYL  (14)  )  ) 


=  4.34  29  4  48*CLOG(CTL  (14)  ) 


=4.34  29  4 48 *C LOG (ZXL  (14) ) 


=  4.34  294  48*3 LOG  (Z 7  L  (14)  ) 


=ATAN2(AIMAG  (COL  (I  4)  )  ,  REAL  (CDL  (14)  )  ) 


COL  (14)  =CSQRT(COL  (14)  *CONJG  (COL  (14)  )  ) 


COPOL(I4)=SPOL(I4)/  (CS2RT  (SPO  (14)  )  *CSQRT  (SPL(I4)  )  ) 


COPOL(I4)  =CSQRT  (COPOL(I4|  *CONJG  (COP3L  (14)  )  ) 


COHOI(IU)  =SMOL(I4)/(CSQRT  (SMO  (14)  )  *CSQRT  (SML  (14)  )  ) 


C0M0L(I4)  =CSQRT  (COMOL  (14)  *CONJG  (COMDL  (14)  )  ) 
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SPO  (14)  =4.34  294  48*CL0G(SP0  (14)  ) 

SPL  (14)  =4.3429!*48*CLOG(SPL  (14)  ) 

SMO  (14)  =4.34294 48*CLOG(SMD  (14)  ) 

SML  (14)  =  4.34  29'4  48*CLOG(SML  (14)  ) 
44       CONTINUE 

VER3ATEC    PLOT    OF    CALCULATED    QUANTITIES 

NPTS=10./DELTAF+1. 

•NPTS1    DETERMINES    NUMBER    OF    POINTS 
NECESSARY    TO    PLOT    THE    0    TO     10    HERTZ    RANGE. 

ITB  (2)  =0 

ITB  (3)  =2  0 

ITB  (4)  =10 

I TB  (  6)  =  1 

ITB  (7)  =  1 

ITB  (  12)  =1 

RTB  (  1)=0 
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RTB  (2)=0 


RTB  (3)=ALAB(1> 


READ  (5r3000)  TITLE 


CALL    DRAWP(N?rS,FRQ2,CU03,ITB,RTB) 


ITB  (6)  =1 


RTB  (3)=ALAB(2) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS,FRQ2,CV03,ITB,RTB) 


ITB  (6)  =  1 


RTB  (3)  =ALAB(3| 


READ  (5  r  3  000)  TITLE 


CALL    DRAWP(N?rS , FREQ, CPCO  ,  ITB, RTB) 


I TB  (  6)  =  1 


RTB  (3)  =ALA3  (U) 


READ  (5,  3000)  TITLE 


CALL    DRAWP(NPTS ,FREQ  ,  CRCO  ,  ITB,  RTB) 
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ITB  (6)=1 

RTB  (3)  =  A  LAB (5) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPTS ,FREQ,CROO (2)  , ITB, RIB) 

ITB  (6)  =  1 


RTB  (3)  =ALAB(6) 


READ  (5,3000)  TITLE 


CALL    DRAWP(N?TS , FRQ2 , CIOD , ITB, RTB) 


ITB  (6)  =1 


RTB  (3)  =ALAB(7) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS ,FREQ,CWO0 ,IT3, RTB) 


ITB  (6)  =1 


RTB  (3)=ALA3(8) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS , FREQ, COOD , ITB, RTB) 
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I TB  (  6)  =  1 


RTB  (3)=ALAB(9| 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPrS,FRQ2,CULL,ITB, PTB) 


ITB  (6)  =1 


RTB  (3)  =ALAB  (10) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS,FRQ2,CVLL, ITB, RTB) 


ITB  (6)  =  1 


RTB  (3)=ALAB(  11) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPT5 , FREQ, CPLL , ITB, RTB) 


ITB  (6)  =1 


RTB  (3)  =  ALAB(12) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS  r  FREQ, CRLL  , ITB,  RTB) 
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ITB  (6)=1 

RTB  (3)=ALAB(13) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPrS  rFREQ,CRLL  (2)  , ITB, RIB) 

ITB  (6)  =  1 

RTB  (3)=ALAB(14) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPrS,FRQ2,CILL,ITB,RTB) 

ITB  (6)  =1 

RTB  (3)=ALAB(15) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPTS  ,FREQ,CWLL  ,  I  TB,  RTB) 

ITB(6)  =  1 

RTB  (3)=ALAB(16) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPrS , FREQ, COLL , ITB, RTB) 
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ITB  (6)  =1 


RTB  (3)  =  ALAB(17> 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPrS , FRQ2 , SPOD , ITB, RTB) 


ITB  (6)  =1 


RTB  (3)  =ALAB(13) 


READ  (5,3000)  TITLE 


CALL    DRAWP(N?IS  ,FRQ2  ,  SPLL  ,ITB,  RTB) 


ITB  (6)  =  1 


RTB  (3)  =ALAB(1Q) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPrS , ERQ2 , SMOD , ITB, RTB) 


I TB  (  6)  =  1 


RTB  (3)=ALAB  (20) 


READ  (5,  3000)  TITLE 


CALL    DRAWP(NPTS ,FRQ2 ,SMLL , ITB, RTB) 
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ITB  (6)=1 

RTB  (3)  =  A  LAB (21) 
READ  (5,3000)  TITLE 

CALL    DRAWP(N?TS  ,  FREQ,COPOLL  ,IT  E,  RTB) 
ITB  (6)  =1 
RTB  (3)=ALAB(22) 
READ  (5,3000)  TITLE 

CALL    DRAWP(NPIS , FREQ, COMDLL , ITB, RTB) 
3000      FORMAT(6A8) 
STOP 
END 

/* 

//GO. SYS  IN    DD    * 

TEST    POWER     SPECTRAL    DENSITY    OF    SIN    UHZ 
IN    XXO,30    AVGS,(REF    DB    VS    LOG    FREQ) 
TEST    POWER     SPECTRAL    DENSITY    OF    SOS    UHZ 
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IN  YYO, 30  AVGS,(REF  DB  VS  LOG  FREQ) 
TEST  STOKES  1/STOKES  3  OF  SIN  UHZ 
IN  XXO  AND  COS  UHZ  IN  YYO, 30  AV3S 
TEST  STOKES  2/STOKES  3  OF  SIN  UHZ 
IN  XXO  AND  COS  UHZ  IN  YYO, 30  AV3S 
TEST  STOKES  3/STOKES  0  OF  SIN  UHZ 
IN  XXO  AND  COS  UHZ  IN  YYO, 30  AV3S 
TEST  STOKES  0  OF  SIN  UHZ  IN  XXO  AND 
COS  UHZ  IN  YYO, 30  AV3S,(REF  DB  7S  LOG) 
TEST  PHASE  OF  SIN  UHZ  IN  XXO 
AND  COS  UHZ  IN  YYO, 3D  AVG5 
TEST  COHER  OF  SIN  UHZ  IN  XXO 
AND  COS  UHZ  IN  YYO, 30  AVGS 
TEST  POWER  SPECTRAL  DENSITY  OF  SIN  UHZ 
IN  XXL, 30  AVGS,(REF  DB  VS  LOG  FREQ) 
TEST  POWER  SPECTRAL  DENSITY  OF  COS  UHZ 
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IN  YYL,30  AVGS. (REF  DB  VS  LOG  FREQ) 

TEST  STOKES  1/STOKES  0  OF  SIN  UHZ 

IN  XXL  AND  COS  4HZ  IN  YYL, 30  AV3S 

TEST  STOKES  2/STOKSS  3  OF  SIN  4HZ 

IN  XXL  AND  COS  4HZ  IN  YYL, 30  AV3S 

TEST  STOKES  3/STOKES  0  OF  SIN  4HZ 

IN  XXL  AND  COS  4HZ  IN  YYL,  30  AV3S 

TEST  STOKES  0  OF  SIN  4  HZ  IN  XXL  AND 

COS  4HZ  IN  YYL, 30  AV3S  ,  (REF  DB  VS  LOG) 

TEST  PHASE  OF  SIN  U HZ  IN  XXL 

AND  COS  4HZ  IN  YYL, 30  AVGS 

TEST  COHEE  CF  SIN  4HZ  IN  XXL 

AND  COS  4HZ  IN  YYL, 30  AVGS 

TEST  RT  CIRC  POLARIZATION  PSD  OF  SIN  U HZ  IN 

XXO  AND  COS  UHZ  IN  YYD ,  (DB  VS  LOG), 30  AVGS 

TEST  RT  CIRC  POLARIZATION  PSD  OF  SIN  UHZ  IN 
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XXL    AND    COS    4HZ    IN    YYL ,     (DB    VS    LOG), 30    AVGS 
TEST    LEFT    CIRC    POLARIZ ATION    PSD    OF    SIN    4HZ 
IN    XXO    AND    COS    4HZ    IN    YYO     (DB    VS    LOG)  ,30    AVGS 
TEST    LEFT    CIRC    POLARIZATION    PSD    OF    SIN    4HZ 
IN    XXL    AND     COS    4HZ    IN     YYL     (D3    V3     LOG) ,30    AVGS 
TEST    COHER     FT    CIRC    POLARIZATION    OF    SIN    4HZ 
IN    XXO, XXL     AND    COS    4HZ    IN    YYO,YYL,     30    AVGS 
TEST    COHER     LEFT    CIRC    POLARIZATION    OF    SIN    4HZ 
IN    XXO, XXL     AND    COS    4BZ     IN    YYO, YYL,     30    AVGS 
/* 
// 
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APPENDIX    B 

MASS    SI3RAGE    PROGRAM 
//J0EB1MB    JCB     (1592,01 29) , 'FISHER    SMC  .  1399 •  fCLASS  =  F 

//*MAIN    0RG=NPGVM1.  159  2P 

//    EXEC    FORTXCLG 

//FORT.SYSIN   DD    DSN=MS S. SYS3. NONIMSL. SOURCE (FOURT) , DISP= SHR 

//  DD       * 

C  THE    ARRAY    'IN'     WILL    BE    USED    TO 

C  RECEIVE    THE    DATA    PASSED    FROM 

C  THE    SUBROUTINE    'READ1     AND    THEN 

C  TRANSFERRED    TO    THE    APPROPRIATE 

C  XXX    OR    YYY     ARRAY. 

INTEGER*2    IN(15) 

COMPLEX*3    XXX  (81  92),  YYY(8192) 

DIMENSION    FREQ(3  19  2)  , TIME (8 192)  ,  WORK  (16384) 

DATA     XXX,YYY/153 8U*(0.0,0.  D) / 

DATA    TIME, FREQ/1 6384*0.0/ 
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THE    FOLLOWING    SECTION    READS    THE 
FIRST    SIX    SECONDS    OF    COMPUTER    TAPE 
AND    DISCARDS    THIS    DATA. 

ITL=192 

DO    55    JJ=1 ,ITL 

CALL     RD(20,IN,20 0,IREC,IRR) 
55    CONTINUE 

IFRAME=8192 

NR=2  0 

FNR=  FLOAT  (NR) 

DO    70    L1=1,NR 

THE    NEXT    LOOP    READS    THE    COMPUTER 
TAPE    USIN3    THE    PROGRAM    PROVIDED 
BY    MR.    US     STANTON    OF    NAVAL    POST- 
GRADUATE   SCHOOL. 

DO    60    JJ=1,IFRAME 
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CALL     RD(20,IN,10 00,IREC,IRR) 

XXX(JJ)=IN(2) 

YYY(JJ)=IN(3) 
60    CONTINUE 

N=8192 

FN=FLOAT  (N) 

DELTAT=1./32. 

T=  FN*  DEL  TAT 

DELTAF=1./T 

DO    2  0    J=1  ,N 

TIME  (J)=DELTAT*FLOAr  (J) 

FREQ  (J)  =DELTAF*FLOAI  (J) 

THE    NEXT    I*     STEPS    CONVERT    THE 

DATA    TO    VDLTAGE    AND    ENSURES    THAT 

NO    ERRONEOUS    DATA    HAS    BEEN    INTRODUCED 

INTO    THE    ARRAYS. 
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XXX(J)=*  (XXX(J)-2  04  5.5)*10./2045.5 
XXX<J)=REAL(XXX  (J)  ) 

YYY(J)= (YYY(J) -2  04  5.5) *10-  /2045.5 
YYY(J)=REAL(YYY(  J)  ) 

20    CONTINUE 

CALL     FOURT(XXX,N  ,1,-1 ,0, WORK) 

CALL     FOURT  (YYY,N  ,1,-1,0,  WDRK) 

DO    40    K4=1,N 

XXX(K4)  =XXX(K4)  /FN 

YYY(K4)  =YYY(K4)  /FN 
40    CONTINUE 

THE    NEXT    LOOP    APPLIES    THE    SYSTEM 
TRANSFER    FUNCTION    TQ    THE    TRANSFORMED 
FREQUENCY    DOMAIN    DATA.     THE    TRANSFER 
FUNCTION    CONVERTS    VOLTS    TO    NANOIESLAS 

DO    9     L= 1  ,  N 
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FRQ=FFEQ  (L) 


1    IF  (FRQ.LE.15.)  30    TO    2 


XXX (L)=X XX  (L)  /(1 05.5-3.  14*FRQ) 


YYY(L)=YYY(L)  /(1  31  .3  2-7.  58  8*FRQ) 


GO    TO    8 


2    IF  (FHQ.LE.  10.)  30    TO    3 


XXX(L)=XXX  (L) /(5 . 958*FRQ-3  0. 97) 


YYY(L)=YYY  (L)  /  (7  . 1  66  *FR2-3  9  .  99) 


GO    TO    8 


3    IF  (FRQ.LE.7.5)  20    TO    U 


XXX(L)=XXX  (L)  /(3  .U92*FR2-5.31) 


YYY(L)=YYY(L)  /(4 . 252*FRQ- 1 0 .  8  5) 


GO    TO    8 


4    I?  (FRQ.LE.5.)  GO    TO    5 


XXX(L)=XXX  (L)  /(2  .6311*FKQ*0.  146  67) 


YYY(L)=YYY(L)  /(3  .  0  1  2*FRQ-  1 .  5  5) 
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GO    TO    8 


5    IF    (FBQ.LE.3.)  3D    TO    6 


XXX(L)=XXX(L) /{2 .6311*FRQ+0. 146  6  7) 


7    YYY(L)=YYY  (L)  / (2  .702*FRDJ 


GO    TO    8 


6    XXX(L)=XXX(L) /(2  .72*FRQ) 


GO    TO    7 


8    CONTINUE 


9    CONTINUE 


THE    NEXT    WRITE    STATEMENTS    SEND 


THE    CONVERTED    DATA    ID    MSS 


FOR    FUTURE     MANIPULATION    AND    RECALL, 


THIS    NEXT    SET    OF    STATEMENTS    ALSO 


PROVIDES    THE    USER    A    DIAGNOSTIC. 


WRITE  (21)  XXX 


WRITE  (6,*)XXX(1)   ,XXX(8192> 
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WRITE  (21)  YYY 

WRITE  <6,*)YYY(1)   ,YYY(8192» 
70    CONTINUE 
ENDFILE    21 
STOP 


END 


SUBROUTINE   PD(I UN, 10, IRS, IREC,IRQ) 


THIS     PROCEDURE    FURNISHED    BY    MR.    TIM    STANTON, 


DEPARTMENT    OF    OCEANOGRAPHY. 


READ    DATA    FROM    ION,     ALLIGN    ,     CHECK    &    RETO RN 


ION=TAPE    NOMBER,    EG    20 
I0=INTEGER*2    ARRAY,     16    LONG, 
(VALUES    0-4095  ,    SUBTRACT    2048)  *5/2028 .    GIVES    VOLTAGE 
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IRS=    NUMBER    OF    PESINCS    ALLOWED     (ERRORS) 
IREC=    COUNTER     OF    RECORDS     (FRAMES    OF    DATA) 

BLOCK   512     BITS,  32    BITS    =    RECORD 

800    BPI    TAPE   [JNLABLED 
IRQ=   NUMBER    OF     ACTUAL    RESINCS     (ERRORS) 


INTEGER    *    2    ID  (  16)  ,  IP  (16) 

DATA    IRR    /0/ 
IF     (IREC.EQ.O)     IS=D 
IER=0 
20  FORMAT     (16A2) 

IF     (IS.NE.O)     3D    TO    50 
READ     (IUN/20rEND=900)     IP 
IREC=IREC+1 

ao        is=is+i 


ill 


IF    (IS.  LT.17)     GO    TO    50 
READ    (IUN,20,END=900)     IP 
IS  =  1 

IREC=IREC-H 
50         ICH=IMASK  (IP(IS)  ,3,0) +1 
:  WRITE    (6r55)     ICH,IS,  IUN,IREC 

55  FORMAT     ('       RESYNCING    ICH,  IS  ,IUN ,  IRSC     »  ,418) 


IF     (ICH.NE.1)     GO    TO    40 
DO    100   1=1,16 

IO(I)=ISHIFT  (IP  (IS)  ,4) 
ICH* I  MASK  (IP  (IS)  ,3,0)  +1 
IF     (ICH.EQ.I)     30    TO    80 
IER=IER+1 

WRITE     (6,70)        IUN,IREC    ,I,ICH,IER 
70  FORMAT     ( ■ UNIT' , 13, ■ RECORD ', 16 ,' CHAN    &    DATA    CH    ',214, 
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•ERRORS     ',17) 


80  IS=IS+1 


IF     (IS.  LT.17)     20   TO    103 


READ     (IUN,20,  END=900)     IP 


IS=1 


IREC=IRE01 


100         CONTINUE 


IF    (IER.  EQ.9)     30    TO     150 


IRR=IRR+1 


IF     (IRR.LT.IRS)      GO    TO    123 


WRITE    (6,110) 


110  FORMAT     (M    STOPPED    IN    SUB    RD    BECAUSE 


C    OF    IRR.GT.1  ,15,  •    kl    L113') 


IRQ=IRR 


STOP 
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120         CONTINUE 

WRITE    (6,130)     IREC,  IRR 
130         FORMAT     ('RESYNC    AT    FRAME' , 16 , ' WITH    TOTAL    ERRORS', 17) 

IER=0 

IRQ  =  IRR 

GO    TO   50 
150         CONTINUE 

RETURN 
900         WRITE    (6,910)     I  UN,I REC 
910         FORMAT     ('1         END    OF    UNIT    ',13,'     AT    REC    ',17) 

STOP 

END 

FUNCTION    ISHIFT     (IN,NPLC) 
:  RETURNS    SHIFTED    VALUE    OF    1*2    WORD    IN 

:  -VE    LEFT,  +  VE    RIGHT    SHIFT 
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INTEGER    *    2    IN 


IP=IN 


IF     (IP.LT.O)     IP=IP+65536 

IF    (NPLC.LT.  0)      GO    TO    30 

ISHIFT=IP/(2**I  ABS(NPLC)  | 

RETUBN 
30  ISHIFT=IP*(2**I ABS(NPLC)  ) 

IF    (ISHIFT.GT. 65535)     ISHIFT  =  MOD  (ISHIFT , 65536) 

RETURN 

END 

FUNCTION    IMASK     (IN,  IBL,  IBR) 
:  MASK    1*2    WORD    IN    OUTSIDE    BITS    IBL    &    IBR 


INTEGER    *    2    IN,  10 


10=  IN 


IF     (IBR.FQ.O)     GO    T3    50 
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IT=ISHIFT  (IN,IBR) 
10=1  T 
50  IP=ISHIFT(I0,IBL-15-IBR) 

IO=IP 

IMASK=ISHIFT  (13  ,  15-IEL) 

RETURN 

END 

/* 

//GO.FT21F001  DD   UNII  = 3330 V, MSV3P=PUBUAr DISP= (NEW , CATLG)   , 

//  DSN=MSS.  S1592.GMDI1 A, 

//  DCB=(RECFM=VBS,BLSSIZE=4096, LRSCL=4092)  , 

//  SPACER  YL,(8,8>  ) 

//GO.FT20F001  DD    UN  11=  3400  -4  ,  VOL  =  S  ER=GMDT1  Ar  DISP=  (OLD,  PASS)  , 

//  LA3EL=(1  ,NL,  ,IN)  , 

//  DCB=(RS3FM=FB,LRECL=3  2,BLKSIZE=512,DEN=2) 

// 
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APPENDIX    C 

N  MAIN    PROGRAM 

//JOEK1MB   JOB     (1592,01 29) , 'FISHER    SMC. 1 399 « , CLASS= F 

//*MAIN    0RG=NPGVM1.  159  2P,LINES= ( 60) 

//♦FORMAT    PR,DDNAME=PLOT.SYSVECrR,DEST= LOCAL 

//    EXEC    FRTXCLGP,PARM.  LKED=«LISr, M AP, XREF' , REGION. GO=270 OK 

//SYSUT1  DD       UNIT=SYSDA,  SPACE=  (CYL,  (8,8)  ) 

//SYSUT2         DD       UNIT=SYSDA,  SPACE  =  (CYL,  (8,8)  ) 

//SYSLIN  DD       SPACE= (6080 ,  (80, 30)  )  ,UNIT  =  SYSDA 

//FORT.SYSIN   DD    DSN=3S S. SYS3. NONIMSL . SOURCE (FOURT) ,DISP= SHR 

//  DD       * 

C  DETAILED    DESCRIPTIONS    OF    THE 

C  VARIOUS    STEPS    AND    LOOPS    ARE 

C  CONTAINED    IN    APPENDICES    A    AND    B. 

INTEGER*2    IN  (15) 

COMPLEX*8    XXO(8  192)  ,YYO  (3  192)  , 
C    CSO  (8192)  ,COO(3  192)  ,  ZX0  (8192)  , 
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C    CRO  (819  2)  ,CTO(8  19  2)  ,CUO  (3  19  2)  , 


C    CVO  (8192)  ,CWO(8  192}  ,CP3(3192)  , 


C    ZYO  (8192)  ,COMOL  (819  2)  ,COL  (8  19  2)  , 


C    CUL  (819  2)  rCVL(8  192)  ,CWL  (319  2)  , 


C    COPOL(8192)  ,  XXX  (819  2)  ,YYY  (8  192)  , 


C    XXL  (8192)  ,YYL(3  192)  ,ZXL(3192)  , 


C    ZYL  (8192)  ,CSL(3  192)   ,BP3(8192), 


C    BMO  (8192)  ,BPL(3  192)  ,BML(3192)  , 


C    SPOL  (8192), SMDL  (8192)  ,SPD  (8192)  , 


C    SPL  (819  2)  rSHO(8  19  2)  ,  CTL  (8  19  2)  , 


C    CPL  (8192)  ,CRL(3  192)  ,SML(3192) 


DIMENSION   TIME(8192)  ,  FRE2  (8192)  , 


C    WORK  (16384)  ,  FR2  2(8192) 


INTEGER*4    ITB(12)/12*0/ 


REAL*4  RIB (2  8)  /2  8*0.  0/,  CROO  (16  384)  , 


C    CPOO  (16384)  ,  CU30(16384|  , 
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C    COOO  (16  3  84)  ,  CT3  0(16  3  84)  ,:V00  (16  384)  , 

C    CWOO  (16384)  rCrLL(16384|  , 

C    CPLL  (16384)  ,CRLL(1S384)  ,30LL(16384)  , 

C    CULL  (16384)  ,CVLL(1S384)  , 

C    CWLL  (16384)  ,  CDMOLL  ( 16384}  r  COPOLL  ( 16384)  , 

C    SPOO  (16384)  ,SPLL(1S384)  , 

C    SMOO  (16384)  ,SMLL(16384) 

REAL    ALAB  (22) /'  PDXO»  , 'PDZO1  , 'ST  10' , 
C    «ST20' , • ST301,1  STOO1 ^'PH&O' , 
C    ' GAMG« ,  »PDXL',  '  PDYL» ,'SI1L,,'ST2L', 
C    'ST3L»  ,  •  STOL','  PHAL'  ^SAHL1  , 
C    »  SPLO1 , • SPLL»,  f  SMOO1 , 'SMLL'  , •CPOL' , 'CMOL1/ 

REAL*8    TITLE  (12) 

EQUIVALENCE  (TITLE  (1)  ,RTB(5)  )  ,  (CR00(1)  > 
C    CRO  (  1)  )  ,  (CP0O(1)  rCPO  (1}  )  , 
C     (C000(1)  ,C03  (1)  )  ,  (CT00(1)  ,  CTO  (1)  )  ,  (CU00(1)  ,  CUO  (1)  )  , 
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C  (CW00(1)  ,CWO  (1»  )  ,  (CVO0(1>  ,CV0(1)  )  ,  (~TLL(1)  ,CTL(1)  )  , 
C  (CPLL(1)  ,CPL  (1|  )  ,  (CELL  (11  ,  CRL(1)  )  ,  (C0LL(1)  ,  COL  (1)  )  r 
C     (CULL(1)  ,CUL  (1)  )  ,  (CVLL(1)  ,  CVL  (1)  )  ,  (CWLL(  1)  ,CWL(1)  )  , 

C    (CO?CLL(1)  ,C0PDL(1>  )  ,  (COSOLL(1)  , 

C    COMOL(1)  )  ,  (SPDD  (1)  ,SPO(1)  )  , 

C     (SPLL(1)  ,SPL  (1)  )  ,  (SMOO  (1)  ,  SMO  (1))  ,  (SMLL(1)  ,  SML(1|  ) 

DATA    XXO,YYO,CSO,COO,ZXO,  ZYO,CRO, 
C    CTO, CUO, CVO, CW3 ,XXX,YYY, 
C    XXL, YYL,ZXL,ZYL  ,CSL,BPO, B MO , BPL, 
C    BML, SPOL,SMOL,  SPO,SPL,SMD, 
C    SML, CTL,CRL,COL ,CVL, COPOL , CWL, COPOL, 
C    COMOL/294912*(0 .0,0. 0)  / 

DATA    TIM3,FRE3, FRQ2/2U575 *0.  0/ 

ITL=192 

DO    55   JJ=1,ITL 

CALL    RD  (20,IN,200,IRZC,IRR) 
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55      CONTINUE 

IFRAME  =  8192 

NR=20 

FNR  =  FLOAT(NR) 

THE    FIRST    STATEMENTS    OF    LOOP    70 
RECALLS    DATA    PREVIOUSLY    STORED    IN 
THE    IBM    3033    MASS    STORAGE    SYSTEM. 

DO    70   L1=1,NR 

READ  (21)  XXX 

READ  (21)  YYY 

DO    43   II=1,IFRAME 

XXL  (II)  =XXX(II) 

YYL  (II)  =  YYY  (II) 
43      CONTINUE 

DO    6  0  JJ=1,IFRfcME 

CALL    RD  (20, IN, 1 000, IREC,IRR) 
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XXO  (JJ)  =IN(2) 
YYO  (  JJ)  =IN(3) 

60      CONTINUE 
N=8  1  92 
FN=  FLOAT  (N) 
DELTAT=1./32. 
T  =  FN*DELTAT 
DELTAF=1  ./T 
DO    20   J=1,N 

TIME  (J)  =DELT  Ar*  FLOAT  (J| 
FREQ  (J)  =DELTAF*  FLOAT  (J) 
FRQ2  (J)  =ALOG10  (FREQ  (J)  ) 
XXO  (J)  =  (XXO  (J)  -2045.5)  *  5.  /2  04  5.  5 
XXO  (J)  =REAL(XXO  (J)) 
YYO  (J)=  (YYO  (J)  -2045.5)  *5.  /2  04  5.5 
YYO  (J)  =REAL(YYO  ( J)  | 
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20      CONTINUE 

CALL    FOURT(XXO, N, 1, -1,0, WORK) 

CALL    FOURT(IY0, N, 1, -1,3, WORK) 

DO    40   K4=1,N 

XXO  (K4)  =XXO(K4)  /FN 

YYO  (K4)  =  YYO(K4>  /FN 
40      CONTINUE 

DO    9    L=1  ,N 

FRQ  =  FREQ  (L) 

1  IF(FRQ.  LE.15.)  30    TD    2 

XXO  (L)  =XXO  (L)/(-81.36    +17.066*FRQ    -    0 . 5389*FRQ*FR Q) 
YYO(L)  =YYO(L)  /(-6  3.  435    +1  1 . 71 06*FRQ- 0 .  2  3  279*FRQ*F  RQ) 
GO    TO   8 

2  IF(FRQ.  LE.10.)  30    TD    3 

XXO  (L)  =XXO(L)  /(-18.486    +    4.9029*FRQ) 

YYO  (L)  =YYO(L) /(-63.  435    +11.7106*FRQ    -0 .  23279*FRQ* FRQ) 
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GO    TO   8 


3      IF(FRQ.  LE.7.5)  3  0   TO    4 


XXO  (L)=XXO(L)  /(3.94  7*FRQ-9.3  68) 


YYO  (L)  =YYO(L)  /(  4.29  5*FRQ-11  .8  37) 


GO    TO    8 


U      IF(FHQ.LE.5.)  33    TO    5 


XXO  (L)  =XXO(L)  /(2.8105*FP2-0.5  6  57) 


YYO  (L)  =YYO(L)  /(2.80  0*FRQ-.5  0) 


GO    TO   8 


5      XXO  (L)  =XXO(L)  /(2.723*FRQ) 


YYO  (L)=YYO  (L)  /(2.71  7*FRQ) 


GO    TO   8 


8      CONTINUE 


9       CONTINUE 


DO    30  11  =  1, N 


ZXO  (II)  =ZXO(II)  +  (XXO  (II)  *30NJG  (XXO  (II)  )  ) 
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ZYO  (II)  =ZYO(II)  ♦  (YYO  (II)  *CONJG  (YYO  (II)  )  ) 


CSO  (II)  =CSO(II)  +(XXO  (II)  *C0NJG  (YYO  (II)  )  ) 


ZXL  (II)  =ZXL(II)  +  (XXL(II)  *CDNJG  (XXL  (II)  )  ) 


ZYL  (II)  =ZYL(II)  +  (YYL(II)  *CONJG  (YYL(II)  )  ) 


CSL  (II)  =CSL(II|  «•  (XXL (II)  *C0NJG  (YYL(II)  )  ) 


BPO  (II)  =XXO(II)  +  (  (0.  ,1.)  *YYO  (II)  ) 


BMO  (II)  =XXO  (II)  -(  (0.  ,1.)  *YYO  (II)  ) 


BPL  (II)  =  XXL(II)  +(  (0.  ,  1.)  *YYL(II)  ) 


BML  (II)  =  XXL  (II)  -(  (0.  ,1.)  *YYL  (II)  ) 


SPOL  (II)  =SPOL(I  I)  ♦  (BPO  (II)  *CONJG  (BPL  (II)  )  ) 


SMOL  (II)  =SMOL(II)  + (BMO (II)  *CONJG(BML  (II)  )  ) 


SPO  (II)  =SPO(II)  ♦  (BPO  (II)  *CONJG  (BPO  (II)  )  ) 


SPL  (II)  =SPL  (II)  ♦  (BPL(II)  *CONJG  (BPL(II)  )  ) 


SMO(II)  =  SMO(II)  +  (BMO  (II)  *CONJG  (BMO  (II)  )  ) 


SML  (II)  =SHL(II)  ♦  (BML  (II)  *CONJG  (BML  (II)  )  ) 


3  0      CONTINUE 
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70      CONTINUE 

DO  33  13  =  1, N 
ZXO  (13)  =  ZX0(I3)  *T/FNR 
ZYO  (13)  =ZY0(I3)  *T/FNR 
CSO  (13)  =  CS0(I3)  *T/FNR 
ZXL  (13)  =  ZXL(I3)  *T/FNR 
ZYL  (13)  =ZYL(I3>  *T/FNR 
CSL  (13)  =CSL(I3|  *T/FNR 
SPOL  (13)  =SPOL(I  3)  *T/FNR 
SPO  (13)  =SP0(I3|  *T/FNR 
SPL  (13)  =SPL(I3>  *T/FNR 
SMO  (13)  *SM0(I3)  *T/FNR 
SML  (13)  =SML(I3)  *T/FNR 
SMOL  (13)  =SMOL(I  3)  *T/FNR 

3  3      CONTINUE 

DO    4U   14=1, N 
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CTO  (iu; 


=  (ZXO  (14)  +  ZYO  (14)  I  *2./T 


cpo  (14; 


=  (ZXO  (I '4  )-ZYO(I4)  )  *2./(T*CTO  (14)  ) 


CRO  (14' 


=  (4*CS0  (14))  /(CIO  (14)  *T) 


coo  (14; 


=  CSO(I4)  /(CSQRT(ZXD  (14)  )  *CSQRT  (ZYO  (14)  )  ) 


CWO  (14] 


=ATAN2(AIKAG  (COO  (14)  )  ,  REAL  (COO  (14)  )  ) 


COO  (14] 


=  CSQRT(COO(I4)  *C3SJG  (COO  (14)  )  ) 


CTO  (14; 


=  4.34  29  4  48*3 LOG  (CTO  (14)  ) 


CUO  (14; 


=  4.34  294  4  8*CLOG(ZXO  (14) ) 


CVO  (14] 


=  4.34  29  4  4  8*CLOG(Z!T3  (14)  ) 


CTL  (14; 


=  (ZXL  (I4)+ZYL(I4)  |  *2./T 


CPL  (14; 


=  (ZXL  (14)  -ZYL  (14)  }  *2./(T*CTL(I4)  ) 


CRL  (14] 


=  (4*CSL(I4)  )  /(CTL  (14)  *T) 


COL  (14] 


=CSL(I4)  /(CSQRT(ZXL(I4)  )  *CSQRT  (ZYL  (14)  )  ) 


CTL  (14] 


=  4.34  29  4 4 8*C LOG (CIL  (14) ) 


CUL  (  14] 


=  4.3429448*CLOG(ZXL  (14)  ) 


CVL  (14] 


=  4.34  294  48*CLOG(ZYL  (14)  ) 
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CWL  (14)  =  ATAN2(AIKAG  (COL  (14)  )  ,REAL  (COL  (14)  )  ) 
COL  (14)  =CSQRT(COL  (14)  *C0NJG  (COL  (14)  )  ) 
COPOI(IU)  =SPOL(I4)/(CSQRI  (SPO(I4)  )  *CSQRT  (SPL(I4)  )  ) 
COPOL(I4) =CSQRT  (COPOL (14)  *CONJG  (COPOL  (14)  )  ) 
COMOI(I4)=SMOL(IH)/(CS3Rr  (SMO  (14)  )  *CSQRT  (SML(I4)  )  ) 
COMOL(I4) =CSQRT  (COMOL(I4)  *CONJG  (COMOL (14)  ) ) 
SPO  (14)  =4.3429448*CLOG(SPD  (14)  ) 
SPL  (14)  =4.34  294  48*C LOG  (SPL  (14)  ) 
SMO  (14)  =4.3429448*CLOG(SMO  (14)  ) 
SML  (14)  =4.3429448*CLOG(S3L  (14)  ) 
44       CONTINUE 

NPTS=10./DELTAF+1. 
ITB  (2)  =0 
ITB  (3)  =2  0 
ITB  (4)  =10 
ITB  (6)=1 
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I TB  (  7)  =  1 

ITB  (  12)  =1 

RTB  (1)  =0 

RTB  (2)  =0 

RTB  (3)  =ALAB(1) 

READ  (5,  3000)  TITLE 

CALL    DRAWP(NPTS ,FRQ2,CU03 ,ITB, RTB) 

ITB  (6)  =1 

RTB  (3)  =ALAB(2) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPTS,FRQ2,CVCO,ITB,RTB) 

ITB  (6)  =1 

RTB  (3)=ALAB(3) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPTS , FREQ,CPOD ,ITB, RTB) 

ITB  (6)=1 
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RTB  (3)  =  ALAB(4) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPIS , FREQf CRO0 , I TB, RTB) 


ITB  (6)  =1 


RTB  (3)  =ALAB(5) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPIS ,FREQ,CRCO  (2)  , ITB, RTB) 


ITB  (6)  =1 


RTB  (3)  =ALAB(6) 


READ  (5,  3000)  TITLE 


CALL    DRAWP(NPIS ,  FRQ2  ,CTCO  ,  ITB,  RTB) 


ITB  (6)  =  1 


RTB  (3)=ALAB(7) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS , FREQ, CH CO , I TB, RTB) 


ITB  (6)  =1 
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RTB  (3)=ALAB(8) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS ,FREQ ,C0CO ,ITB, RTB) 


ITB  (6)=1 


RTB  (3)  =ALAB(9) 


READ  (5,  3000)  TITLE 


CALL    DRAWP(NPTS , FRQ2 ,CULL ,ITB, RTB) 


ITB  (6)=1 


RTB  (3)=ALAB(10) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPIS,FRQ2,CVLL,ITB,RTB) 


ITB  (6)  =1 


RTB  (3)  =A LAB  (11) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS  ,  PRE  Q,CPLL  ,1TB,  RTB) 


ITB  (6)  =  1 
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RTB  (3)  =  ALAB(12) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPrS , FREQ, CRLL , ITB, RTB) 


ITB  (6)  =1 


RTB  (3)=ALA3(13) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPrS ,FREQ,CRLL (2)  , ITB, RTB) 


ITB  (6)  =1 


RTB  (3)  =ALAB  (1^) 


READ  (5,  3000)  TITLE 


CALL    DRAWP(NPIS ,FRQ2 , CTLL , ITB, RTB) 


ITB  (6)  =  1 


RTB  (3)  =  ALAB(15| 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS , FREQ, CWLL , I TB, RTB) 


ITB  (6)  =1 
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RTB  (3)=ALAB(16) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPTS  , PRE Q,COLL ,ITB, RTB) 


ITB  (6)  =1 


RTB  (3)  =ALAB(17) 


READ  (5r  3000)  TITLE 


CALL    DRAWP(NPIS  , FRQ2  ,  SPCO  ,  ITB ,  RTB) 


ITB  (6)  =  1 


RTB  (3)=ALAB(18) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPrS , FRQ2 , SPLL  ,  ITB, RTB) 


I TB  (  6)  =  1 


RTB  (3)  =  ALAB(19) 


READ  (5,3000)  TITLE 


CALL    DRAWP(NPrS ,FRQ2 ,SM0D , ITB, RTB) 


ITB  (6)  =1 
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RTB  (3)  =ALAB(20| 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPrS , FRQ2 , SHLL , ITB, RTB) 

ITB  (6)  =1 

RTB  (3)=ALAB(21) 

READ  (5,3  000)  TITLE 

CALL    DRAWP(NPTS  , FREQ, C0P3LL , ITB, RTB) 

ITB  (6)=1 

RTB  (3)  =ALAB(22) 

READ  (5,3000)  TITLE 

CALL    DRAWP(NPTS ,FREQ ,COMOLL , IT B, RTB) 


3000       FORMAT  (6  A8) 


STOP 


END 


SUBROUTINE    RD (I  UN, 10 , IRS, IREC,IRQ) 
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THIS     PROCEDURE    FURNISHED    BY    MR.    TIM    STANTON, 


DEPARTMENT    OF    OCEANOGRAPHY. 


READ    DATA    FROM    UN,     ALLIGN    ,     CHECK    &    RETURN 


IUN=TAPE    NUMBER,    EG    23 
I0=INTEGER*2    ARRAY,     16    LONG, 

(VALUES    0-4095  ,    SUBTRACT    2048) *5/2028 .    GIVES    VOLTAGE 
IRS=   NUMBER    OF     RESINCS    ALLOWED     (ERRORS) 
IREC=    COUNTER     OF    RECORDS     (FRAMES    OF    DATA) 

BLOCK    5  12     BITS,         32    BITS    =    RECORD 

800    BPI    TAPE    UNLABLED 
IRQ=   NUMBER    OF     ACTUAL    RESINCS     (ERRORS) 


INTEGER    *    2    10  (  16)  ,  IP  (1  6> 
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DATA    IBB    /O/ 


IF     (IBEC.EQ.O)     IS=0 


IER=0 


20  FORMAT    (16A2) 


IF     (IS.NE.O)     3D    TO    50 


READ     (IUN,20,END=900)     IP 


IREOIBE01 


40  IS=IS+1 


IF     (IS.  LT.17)     30    TO    50 


READ     (IUN,20, SND=900)     IP 


IS  =  1 


IREC  =  IREC  +  1 


50  ICH=IMASK  (IP(I3)  ,3,0)  +1 


WRITE    (6,55)     I3H,IS  ,IUN,IREC 


55  FORMAT     ('       RESYNCING    ICH,  IS  ,  IUN ,  IRE:     ',418) 
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IF     (ICH.NE.1)     33    TO    UO 


DO    100  1  =  1,16 


IO(I)  =  ISHIFT  (IP  (IS)  ,4) 


ICH=IMASK  (IP(IS)   ,3,0) +1 


IF    (ICH.  EQ.I)     S3    TO    80 


IER=IER+1 


WRITE     (6,70)        I[JN,IREC    ,I,ICH,IER 


70  FORMAT     (• UNIT' , 13 ,' RECORD « , 16 ,« CHAN    S    DATA    CH    ',214, 


•ERRCRS     »,I7) 


80  IS=IS+1 


IF     (IS. LT.17)     30    TO    103 


READ     (IUN,20,END=900)     IP 


IS=1 


IREC=IREC  +  1 


100         CONTINUE 
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IF     (I3R.EQ.0)     30    TO     150 

IRR  =  IRR+1 

IF     (IRR.LT.IRS)     GO    TO    120 

WRITE    (6,110) 
110         FORMAT     (M    STOPPED    IN    SUB    RD    BECAUSE 
C    OF    IER.GT.»,I6f  ■    AT    L110«) 

IRQ=IRR 

STOP 
120         CONTINUE 

WRITE    (6,130)     IRECfIRR 
130  FORMAT     ('RESYNC     AT    FRAME f , 16 ,  • WITH    TOTAL    ERRORS', 17) 

IER=0 

IRQ  =  IRR 

GO    TO   50 
150         CONTINUE 

RETURN 
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900         WRITE    (6,910)     IUN,IREC 


910         FORMAT     (M         END   OF    UNIT    ',13,'     AT    REC    »rl7) 


STOP 


END 


FUNCTION    ISHIFT     (IN,NFLC) 


RETURNS    SHIFTED    VALUE    OF    1*2    WORD    IN 


-VE    LEFT,+VE    RI3HT    SHIFT 


INTEGER    *    2    IN 


IP=IN 


IF     (IP.LT.O)     IP=IP+65536 


IF     (NPLC.LT.  0)      GO    TO    30 


ISHIFT=IP/(2**I  ABS(NPLC)  I 


RETURN 


30  ISHIFT=IP*(2**E  ABS(NPLC)  ) 


IF     (ISHIFT. GT. 65535)     ISHI  FT  =  MOD  (ISHIFT  ,  65536) 
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/* 


RETURN 


END 


FUNCTION    MASK     (IN,IBLrIBR) 


MASK    1*2    WORD    IN    OUTSIDE    BITS    IBL    &    IBR 


INTEGER    *    2    IN,  10 
IO=IN 

IF     (IBR.  EQ.O)     30    TD     50 
IT=ISHIFT(IN,IBR) 
10=  IT 
50  IP=ISHIFT(I0,IBL-15-IBR) 

10=  IP 

IMASK=ISHIFT  (13  ,15-IBL) 
RETURN 
END 
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//GO.FT20F001    DD   UNIT  = 3400-4 ,  V3L  =  SER=GMDT1 B, DISP= (OLD, PASS) , 

//  LABEL=(1  ,NL,  rIN)  , 

//  DCB=(RECFM=FB,LRE:L=3  2,BLKSIZE=512,DEN=2) 

//GO.FT21F001     DD   UN  11= 3330 V, MSV3P=PUB4 A, DISP= (OLD , KEEP)  , 

//  DSN=MS5. S1592.GMDT1A, 

//  DCB=(F.E3FM=VBS,BLKSIZE=4  096,LRECL  =  40  92) 

//GO.SYSDOMP  DD    SYSOUI-k 

//GO. SYS IN     DD    * 

X    COIL    PSD,     MIRY    BAY, 1  122- 1254    L3CAL , 17AUS82 , 

20    AVGS,    32S/SEC,5VOLT  ,     (DB    NT**2    VS    LOG    FREQ) 

Y    COIL    PSD,    MTRY    BA Y,  1  122- 1254    LOCAL, 1 7AOG82, 

20    AVGS,    32S/SEC,5VOLT ,     (DB    NT**2    VS    LOG    FREQ) 

STOKES    NO.     ONE,     MTRY    BAY,     1122-1254 

LOCAL,     17AUG82,    20    AV3S,    3  2S/SE3,     5VOLT 

STOKES    NO.     TWO,     MTRY    BAY,     1122-1254 

LOCAL,     17AUG82,    20    AV3 S,    3  2S/SEC,     5VOLT 
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STOKES    NO.     THREE,    MTRY    BAY,     1122-1254 
LOCAL,     17AUG82,    20    AVGS,    32S/SEC,    5VOLT 
STOKES    NO.     ZERO,    fURY     BAY,     1122-1254    LOCAL, 
17ADG82,    20    AVGS,    32S/SEC,    5VOLI,     (LOG    VS    LOG) 
PHASE    OF    X     AND    Y,    MTRI    BAY     1122-1254 
LOCAL,     17AUG82,     20    AV3S,    32S/SE3,     5VOLT 
COHER    OF    X     AND    Y   COILS,    MTRY    BAY,     1122-1254 
LOCAL,     17AUG82,     20    AVGS,    3  2S/SE3,     5VOLT 
X    COIL    PSD,    LA    MESA, 11 22-1 2545  LOCAL, 17AUG82, 
20    AVGS,    32 S/SEC,10VOLT,     (DB    NT**2    VS    LOG    FREQ) 
Y   COIL    PSD,    LA    MESA, 1 1 22-1254    LOCAL, 1 7ADG82 , 
20    AVGS,    32S/SEC,10VOLT,     (DB    NT**2    VS    LOG    FREQ) 
STOKES    NO.     ONE,    LA    MESA,     112201254 
LOCAL,    17ADG82,    20    AV3S,    32S/SE3,     10VOLT 
STOKES    NO.     TWO,    LA    MESA,     1122-1254 
LOCAL,     17AUG82,    20    AV3S,    32S/SEC,     10VOLT 


142 


STOKES    NO.     THREE,    LA    MESA,     1122-1254 
LOCAL,     17AUG82,    20    AV3S,    32S/SEC,     10VOLT 
STOKES    NO.     ZERO,    LA    MESA,     1122-125U    LOCAL, 
17ADG82,    20    AVGS,    32S/SEC,     10VOLT,     (LOG    VS    LOG) 
PHASE    OF    X     AND    Y,    LA    MESA,     1122-1254 
LOCAL,     17AUG82,     20    AV3S,    32S/SE:,     10VOLT 
COHER    OF    X     AND    Y   COILS,    LA    MESA,     1122-1254 
LOCAL,     17AOG82,     20    AV3S,    32S/SE3,     10VOLT 
RT    CIRC    POLARIZATION    PSD,MTRY    BA Y , 1 1 22-1 254 
LOCAL, 17AUG82, 20    AVGS, 32S/SEC, 5V3LT, (DB    VS    LOG) 
RT    CIRC    POLARIZATION    PSD,    LA    MESA,     1122-1254 
LOCAL, 17AUG82, 20    AVGS,  32S/SEC, 1 0  VOLT,  (DB    VS    LOG) 
LEFT    CIRC    POLARIZATION    PSD,    MIR?    BAY, 1  122-1 2 54 
LOCAL, 17AUG82,     20    AV3S  ,32S/SEC , 5  VOLT,  (  DB    V5    LOG) 
LEFT    CIRC    PCLARIZAIIDN    PSD,    LA    SESA,     1122-1254 
LOCAL, 17AOG82, 20    AVGS,  32S/SEC, 1 0 VDLT, (DB    VS    LOG) 
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COHER    RT    CIRC    POLARIZATION    HTRY    BAY/LA    MESA, 
1122-1254    LOCAL,    17AU3  82,    20    AV3S,    32S/SEC 
COHER    LEFT    CIRC    POLARIZATION    MIRY    BAY/LA    MESA, 
1122-1254    LOCAL,    17AU382,    20    AV3S,     32S/SEZ 
/* 
// 
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APPENDIX   D 
DATA   CURVES 
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